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Abstract. A family of models for fluctuating loops in a two dimensional random 
background is analyzed. The models are formulated as 0(n) spin models with quenched 
inhomogeneous interactions. Using the replica method, the models are mapped to the 
M — > limits of M-layered 0(n) models coupled each other via <f>\^ primary fields. 
The renormalization group flow is calculated in the vicinity of the decoupled critical 
point, by an epsilon expansion around the Ising point (n = 1), varying n as a continuous 
parameter. The one-loop beta function suggests the existence of a strongly coupled 
phase (0 < n < n*) near the self-avoiding walk point (n = 0) and a line of infrared 
fixed points (n* < n < 1) near the Ising point. For the fixed points, the effective 
central charges are calculated. The scaling dimensions of the energy operator and the 
spin operator are obtained up to two- loop order. The relation to the random-bond 
g-state Potts model is briefly discussed. 



1. Introduction 

A variety of random, or disordered systems are known to have non-trivial universalities. 
These include the universality consisting of the well-known three symmetry classes 
which appear in the weak localization effects to the conductance [53] or the level 
correlations [54] in disordered electronic systems, and the universality in the wave 
function multifractality at the integer quantum Hall transition in which disorder plays an 
essential role [55]; even the universality in the distribution of the zeros of the Riemann 
zeta function is expected to have its origin in some chaotic system [45], and it may 
possibly be closely related to some disordered system. It is amazing that some of these 
systems in appropriate limits are phenomenologically well described by such simple 
ansatz as the one in the random matrix theory. It is, however, by no means certain that 
these phenomena are well understood theoretically. The most fundamental problem to 
be understood is how the universalities in disordered systems emerge from a microscopic 
structure of each system. 

A lot of effort has been made to investigate this direction, and in some cases we 
have a partial answer; we can first formulate a system microscopically, then by using a 
very crude approximation (such as 'dimensional reduction' used in disordered electronic 
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systems [54]), deduce universality directly in some carefully chosen limit. Even under 
such fortunate circumstances, it is usually the case that we do not have much knowledge 
on the way along which the system deviates from the universality away from the limit. 
This is, roughly speaking, because a disordered system formulated on a microscopic 
ground typically has highly nonlinear interactions induced by disorder and becomes 
almost intractable without any crude approximations. Thus, it is important to study 
a simple model which is tractable with a sensible approximation; analysis of such a 
model and ideas used there may lead to a practical way to obtain qualitative results on 
the deviation from the universality in more complicated disordered systems and to an 
insight on the emergence of the universality itself. 

With this broad motivation in mind, among the diverse disordered systems, we 
take a problem in statistical physics, where relations between microscopic properties 
and macroscopic behavior of models are often quite intelligible. More specifically, we 
study a model with quenched disorder defined on a lattice focusing on its possible critical 
behavior. Understanding the behavior of such a system governed by the Hamiltonian 
which contains quenched inhomogeneous interactions is especially important, both from 
the practical perspective that there are no translationally invariant perfect crystals in 
the real systems and from theoretical interests inspired by the other disorder-induced 
phenomena. 

As is well known, the solid notion of the universality class has been established in 
the field of critical phenomena without disorder. Nonlinear interactions inherent in a 
model are essential to understand the existence of the non-trivial universality classes 
[25] ; one should confront with the nonlinearities in a more powerful framework than the 
mean field approximation. The ideas of scaling and renormalization group (RG) give 
such a framework. In this framework, one introduces a theory space and a vector field 
on it, which generates a flow that indicates the scale dependence of a theory considered. 
Assuming critical points in the second order phase transitions are described by scale 
invariant field theories, they correspond to the fixed points of the flow. The flow can 
be determined, step by step, by a path integration on the fluctuations belonging to 
a certain scale; the nonlinear interactions induce the mutual coupling among degrees 
of freedom at different scales. As a result, the global configuration of the flow can 
become non-trivial. Then we can read off the universality classes from the sets of RG 
eigenvalues which characterize the flow linearized around the fixed points. Furthermore, 
it is known that the ideas of RG are not only powerful to resolve the universality classes 
in the critical phenomena of pure systems but also flexible enough to deal with some 
disordered systems. 

For the target of this paper, a system governed by an inhomogeneous Hamiltonian, 
it is natural to consider a configuration of the interactions as one realization from an 
ensemble that respects given probability distribution. Basic observables to be discussed 
can be related to the free energy in a large enough system, which is expected to be 
self-averaging over the disorder distribution. We should therefore study the quenched 
average: the average over the disorder distribution, of the free energy obtained by tracing 
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over the statistical mechanical degrees of freedom. 

Actually, evaluating the quenched average of the free energy is notoriously difficult 
task, since we should evaluate the average of the logarithm of the partition function in 
a non-translational invariant realization of the disorder. It is practical to use either one 
of the two methods: replica [8] or supersymmetry [54, 30]. When these methods are 
applicable, we are left with an effective theory with the translational invariance, but 
this time, with the additional nonlinear interactions induced by the disorder. It should 
be noted that the system acquires enhanced symmetry, namely, the replica permutation 
symmetry or the supergroup symmetry, respectively. 

For a weakly disordered system, one can consider the disorder-induced couplings in 
the effective theory as a perturbation to the corresponding pure theory. A natural 
question is whether the induced coupling destructs the critical phenomena of pure 
system, or not. A basic and general result on this direction is known as the "Harris 
criterion" that tells when the disorder can be neglected at large scales; if the couplings 
are irrelevant in the RG sense, the disorder can not change the universality class of the 
system from that of the pure system [29]. However when they are relevant or marginal, 
one should work harder to analyze the flow, namely, proceed to calculation of loop 
diagrams formed by the disorder- induced couplings. This procedure corresponds to the 
calculation of the RG beta function up to the second or higher order in the couping. 
In the case that the beta function has a zero in addition to the one corresponding to 
the pure fixed point, the flow can (depending on the sign of the beta function) transfer, 
at large scale, the theory to another fixed point; one recognizes that the disordered 
system belongs to a non-trivial universality class. The RG eigenvalues belonging to 
the universality class can be predicted perturbatively. In general, without a special 
reason to be integrable, one can not solve given interacting statistical model. Therefore, 
the importance of such a crossover from one fixed point to the other can not be 
overemphasized . 

We shall study a one-parameter family of disordered models and discuss their 
universality class in two dimensions. Thus, let us briefly note the special role of two 
dimensions in the study of the universality class of pure systems first, and then comment 
on the current status of corresponding study in disordered systems. Working in two 
dimensions provides us with an ideal circumstance to study the crossover in depth. 
First, in two dimensions the conformal symmetry is infinite dimensional, and assuming 
the conformal symmetry (instead of the scale invariance only) enables us to classify 
the possible fixed points of unitary theories with minimal symmetry [3]. Second, it 
can be shown that the RG flow is irreversible for unitary theories ( "c-theorem" ) [22]. 
Third, there are known examples of the crossover in which, at the non-trivial fixed 
point reachable via the flow from other fixed points, the realized enhanced-symmetries 
are known and the non-perturbative results can be obtained [4, 5, 6]. 

In the disordered systems in two dimensions, on the other hand, fixed points created 
by disorder are expected to be described by CFT's, but by more general, non-unitary 
ones. The lack of the unitarity gives rise to the challenging problems of determining 
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the universality class in disordered systems. Major current examples concerned with 
the restrictive cases where the supersymmetry method can be used. They typically 
lead to logarithmic CFT's where the disorder-induced coupling is marginally irrelevant 
[30, 31, 58]. In these cases, the coupling goes to zero under the RG flow, and the models 
do not cause any non-trivial crossover. There is also a consideration on the running of 
the effective central charge defined on the basis of the supersymmetry [31], in analogy 
with the c-theorem in pure (unitary) systems [22]. At present, however, not much is 
known about the crossover cases, and hence our understanding of the disordered critical 
points is quite limited. 

One of a few exceptional examples showing the crossover is the random-bond g-state 
Potts model in two dimension| [8, 9, 20]. The model at q = 2 is the random-bond Ising 
model and does not show the crossover; it is equivalent to the random-mass fermion 
model, or using the replica method, to the multi-color Gross-Neveu model, in all of 
which the interaction is marginally irrelevant [30]. Now in the Z 2 -invariant scalar field 
theory (without disorder), a nontrivial fixed point emerges when the dimension d is 
considered as a continuous parameter in the range d < 4 [25]. In the random-bond 
g-state Potts model, it is also fundamental to consider the parameter q as a continuous 
number. One has then a non-trivial fixed point in the region q > 2. The perturbative 
calculation of the RG eigenvalues belonging to this non-trivial universality class is well 
under control. In particular, the theoretical prediction for the exponent of spin-spin 
correlation function is in good agreement with the numerical simulation [10]. 

In this paper, we introduce the disordered version of the 0(n) model§ and study 
the crossover in it. Our model has its own physical importance at certain integral values 
of n; it corresponds to the random-bond XY (n — 2), the random-bond Ising (n — 1) 
and the polymers (or, self- avoiding walks) in random environment (n — 0). But more 
importantly, we consider the models for continuous values of n; this leads to another 
example of the crossover in disordered system. 

For continuous values of n, as explained below, we have a family of models which 
describe fluctuating loops in a two dimensional random background. As is well known, 
elementary excitations in a spin system like the 0(n) model can be considered as non- 
local geometrical objects, namely, loops in the high-temperature expansion. It has 
been repeatedly emphasized that these loops are analogous to the closed trajectories 
of some particles [43, 44]. Now the parameter n controls a statistics of particles in the 
random potential; we expect distinct, non-trivial behavior according to the value of n. 
In order to study these, we use the replica method as in the studies of the random-bond 
Potts model [8, 9, 20]; the method leads us to consider several, say M, layers of two 
dimensional 0(n) models coupled each other via the disorder-induced coupling, and to 
take the M — > limit in the end. 

The loops in the pure 0(n) model become scale invariant at critical temperature for 
| n | < 2. We investigate critical behavior in the disordered 0(n) model using conformal 

| Other exceptions includes the disordered Dirac fermion problems [59, 60]. 

§ The 0(n) model is another natural extension of the Ising model other than the g-state Potts model. 
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perturbation theory around the one-parameter family of CFT's corresponding to a line 
of the pure 0(n) critical points [12]. The existence of the crossover in our model suggests 
that there is a one-parameter family of non-unitary CFT's. In this respect, we mention 
recent development of the stochastic Loewner evolution (SLE) [32], and its application 
to non-unitary critical points [50, 51]. The line of the fixed points in our model may 
serves as a natural target of such study in disordered system. 

The organization of the paper is the following. In Section 2, we formulate the model 
on a lattice, and explain the types of the quenched disorder considered. Then we use the 
replica method and take the disorder average. As a result, we reach an intriguing picture 
of particles going up and down across the two dimensional layers, thus forming a whole 
connected diagrams. Then we discuss the relation between the observable on a lattice 
and the scaling fields in a continuum limit. We see that the nontrivial cases in which 
the disorder is relevant occur for n < 1. In Section 3, we perform one loop calculation 
and discuss the existence of a non-trivial fixed point. One loop beta function suggests 
the existence of a threshold n*; the fixed point exists for n* < n < 1, while the strongly 
coupled phase exists in < n < n*. In Section 4, we perform the two loop calculation for 
n* < n < 1, using the full information of the four-point function provided by the 0(n) 
CFT. The next-leading order correction for the thermal exponent and the lowest order 
correction for the spin exponent are then found. In Section 5, we calculate the effective 
central charge defined in the replica formalism. We find that this increases, along the 
flow, against the c-theorem which is responsible for unitary theories. We conclude in 
Section 6, and comment on few further directions. Appendix A provides the formulas on 
the critical Liouville field theory used in the paper. In Appendices B and C, we describe 
the calculation of the integral in two loop calculation of the beta function and in the 
spin scaling dimensions, respectively. Finally, Appendix D is devoted to the derivation 
of the integral formula and the expansion techniques. The main formula takes the form 
of the scattering amplitude reflecting the picture that the particle forming the loop can 
propagate via intermediate states while going across the replica layers. 

2. Formulation of the Model 

In this section, we formulate disordered 0{n) loop models on a lattice and discuss 
their continuum limit. In Section 2.1, three types of the disordered lattice models are 
mapped to homogeneous coupled lattice models by the replica method. We introduce 
known CFT for the homogeneous 0(n) model and discuss the relation between lattice 
and continuum in Section 2.2. In Section 2.3, the effective action for the continuum 
disordered model is derived making use of the operator product expansion (OPE). 
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2.1. Disordered Loop and Coupled Loop Model on a Lattice 

We start with the partition function of pure 0{n) model on a two-dimensional lattice: 
Z(t,n)= fY[Ks t )d n Si H(l + t Si - Sj ), (1) 

where Sj is a n-component spin on a site i, and n(s) represents a measure on the 
isotropic internal space; using the notation Tr s . for the tracing operation f Y\ i fi(si)d n Si-, 
it satisfies Tr s . 1 = 1, Tr s . s$ = and Tr s . (sj • Sj) = n. The interaction is short-ranged, 
and the notation refers to a link between the nearest-neighbor sites of the lattice. 

A basic idea in this paper is to continue the parameter n to non-integral values and 
to take advantage of known continuum properties of the corresponding 0(n) model [12]. 
On a lattice, the partition function (1) of the spin system with an arbitrary value of n 
can be interpreted as a model for fluctuating loops. It becomes especially simple on a 
honeycomb lattice, and then tracing over the spin degrees of freedom yields, 

Z(t,n)= Yl t #bonds n #loops , (2) 

config. of loops 

where the summation is taken over the configuration of the closed loops [1]. Now the 
weight per bond of loops is t, while the weight per closed loop is n, which is called a 
fugacity and can take non- integral values. The model in (1) shows universal critical 
behavior when |n| ^ 2, and we may use the same continuum description regardless of 
specific lattice structures ||. 

Another interpretation to make sense out of the non-integral values of n is possible 
for \n\ ^ 2. By orienting the loops, one can assign a complex Boltzmann weight e %x 
(e~* x ) to each clockwise (anti-clockwise) loop with a phase angle x determined from the 
relation n = e %x + e~ %x = 2 cos x- This notion of the oriented loops here is standard in 
the coulomb gas (CG) methods [18], where the loops are mapped to the level lines of a 
height model [18, 15, 33]. One can also consider the loop as a closed trajectory of some 
particle. The local phase factor exp(i^x/27r) is then associated, at a site, with each turn 
of the particle through an angle 9. This is a lattice version of the spin factor discussed 
in [44]. In a sense, the value of n controls statistics of the particles. 

The qualitative behavior of the 0(n) model can be summarized as follows. When 
|n| ^ 2, the model has three different phases separated by a critical point t = t c ^ . 
The region t < t c corresponds to the high temperature phase of the spin model, and 
the length of loop measured by the unit of the lattice spacing is finite. On the other 
hand, the average length of the loop is divergent either at t = t c or in t > t c . The 0{n) 
model at the critical point t — t c is called in the "dilute" phase, since the fraction of the 
number of sites visited by some loop is zero. In t > t c , this fraction becomes non-zero 
and the corresponding phase is called "low temperature" or "dense" phase. It is known 

|| This is valid in the dilute phase but not in the dense phase. For instance, the intersections which 
may occur in the model on a square lattice become relevant in the dense phase and discriminate it from 
the model on the honeycomb lattice. For an explanation of the dilute and the dense phase, see below. 
% On a hexagonal lattice, to be concrete, it is given by t c (n) = (2 + \/2 — n) -1 / 2 . 
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that the shapes of the loops in the dilute and the dense phase are described by the SLE K 
with the Brownian motion amplitude k < 4 and k > 4, respectively [32]. 
We now give the partition function of the random model as 

Z[{t},n] =Tr Si ]J(l + t ijSi . Sj ), (3) 

where the local interactions between spins are position dependent and the notation 
{£} refers to some definite configuration of the interaction. We consider the configuration 
{t} as a realization taken from some ensemble with a probability distribution functional 
P[{t}]. One might assume a non-local probability distribution functional, but here we 
restrict ourselves to study the case of short range correlation. This means that the 
interaction t = tij on each link independently respects single distribution function P(t). 
Later, we shall let P(t) a Gaussian-like distribution. 

Given an inhomogeneous realization {t}, the tracing over the spin degrees of 
freedom is still possible: 

#loops 

Z({t},n)= ^ n* l °° PS Yl *a,(l)a,(2)*a,(2)a,(3) " ' * a, (L, )o,(l), ( 4 ) 

config. of loops 1=1 

where ai(i) (i e {1, 2, • • • , L{\) denotes a site which belongs to the /-th loop of a length 
L[. Now the path of the particles forming loops should avoid the links with higher cost; 
we expect its behavior to change, according to the value of n. The situation is, to some 
extent, analogous to the problem of the electrons in a random potential which has been 
studied in connection with the Anderson localization [52]. As the weight t is different 
from link to link (Figure l-(a)), the exact methods applicable in the pure 0(n) model, 
such as the CG method, can not be applied here. 

Assuming the short-range correlation between the disorder, we can study the self- 
averaging quantities such as the free energy and the translationally averaged correlation 
functions by calculating the quenched average [16]. For example, neglecting the surface 
effects, the free energy of a large enough system A can be considered as a sum of the free 
energies of many macroscopically large subsystems of A each with a different realization. 
Since by the short-range assumption these realizations are independent each other, the 
total free energy per site of A in the thermodynamical limit takes, with the probability 
one, the quenched averaged value defined by 

7 = i im / II dt iJ P (tij) (- In Z [{t}, n]) /N. (5) 

Here, we have used the fact that a free energy per site of a subsystem with N sites 
is given by f — (—\nZ)/N. The overline is used to indicate the averaging over the 
distribution. It should be noted that, in this argument, we need a macroscopic number 
of macroscopically large subsystems. 

In this paper, we use the replica method to evaluate the various quantities related 
to the quenched average of the free energy (5). This method is based on the identity: 

Z M — 1 

\nZ = lim — — — . (6) 
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Using this, the problem of calculating the quenched average In Z is reduced to the task 
of obtaining another quenched average Z M correctly for M ~ 0. To evaluate the latter, 
we first prepare the M G N layers of replicated inhomogeneous 0(n) models with the 
same realization {£}, and then take an average over the distribution. Since, in general, 
the moments of the interaction are non-zero, the resulting theory is a coupled M 
layers of 0(n) models. We calculate the quantity Z M for finite M G N, take the limit 
M — > in the end. This procedure gives, at least formally, gives the desired average. 




Figure 1. (a) Loops on a inhomogeneous lattice. The grayscale on each link represents 
the magnitude of the weight . (b) Loops on several sheets of homogeneous lattices 
are coupled each other. An example of a term with the first order in the coupling ti- 



Although we perform the detailed analysis in continuum theory, let us proceed, for 
a while, within the discrete lattice formulation. The purpose is two-fold: (i) to get some 
insight on underlying physics, as we can concretely see the elementary processes existing 
in the model and (ii) to see that there are possibly many lattice models described by 
the same continuum field theory. Now, we consider the M copies of the model (3) and 
take the disorder average. This gives 

- M 

(z [{ t }, n\) M = j n n u( i+ £ ^ a) • s f) 

M / M M 

= n n 1 + *. e > ■ *t ' + *. e ^ ■ 4°' • •. 

a=l \ a=l a<b 

+ (higher order terms) , 

where in the last line, the first two moments are denoted as t\ = tij and ti = £?■, and 
we omit the terms with the higher moments for simplicity. The term t\ {^s^f 1 ■ 
represents a walk on the same replica layer; the generic case with t\ ^ is discussed in 
the following. In addition to this, the theory acquires nonlinear couplings. For instance, 
when we look at a particular path of a loop, the coupling t 2 ^s\ a ^ ■ s^ sf^ ■ s*pj for 
a < b is considered as a branching process across the replica direction as shown in 
Figure l-(b). On a lattice, evaluating directly the contribution of the diagrams with the 
second or higher orders in such processes should involve some sophisticated enumerative 
combinatorics. Notice, however, that the translational invariance of the individual 0(n) 
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model is now restored. Thus, in the continuum limit, we can use the knowledge of the 
homogeneous 0(n) CFT. 

The restriction on the sum to off-diagonal sector (a < b) in the t<i term follows from 
the fact that the original partition function (3) is linear in the bonds {t^ Sj • Sj). When 
we talk about the universality, however, this formal linearity and resulting off-diagonal 
property should be reconsidered; we should rather consider a larger parameter space 
including more general couplings generated under block spin transformations. 

In order to see this, firstly, it is instructive to remind that the pure model (1) 
corresponds to the Hamiltonian H — — ^2 ln(l + tSi ■ Sj), and originates from the high 
temperature expansion of another spin system which has the same 0(n) symmetry but 
has the different Hamiltonian H = — (3 • Sj. From the RG perspective, these two 
models are regarded as two different points in the same theory space determined from the 
0(n) symmetry. Actually, both systems are believed to belong to the same universality 
class. In fact, with the knowledge of the scaling dimensions of fields obtained by the 
CG method [18], we may argue that these two Hamiltonians are equivalent modulo 
irrelevant operators. In particular, the field (s, • Sj) 2 represents the double occupancy of 
the bonds in the partition function, and the presence of it enables the loops to overlap. 
Although such a field is not contained in the partition function (1), it is natural to think 
that the term (s, • Sj) 2 is generated, under block spin transformation, from the single 
occupation (s, • Sj). But, as we will see in the end of Section 2.2, the term (s, • Sj) 2 is 
irrelevant in the dilute phase, or equivalently, near the critical point. Because of this, 
we can say that the linear expression (1) nicely summarize the characteristics of the 
dilute phase in generic models with the 0(n) symmetry. 

Taking these considerations into account, we now briefly discuss the theory space of 
the coupled models derived from the disordered ones. Our concern is near the decoupled 
critical point; applying the knowledge of the pure 0{n) model, we note the inessential 
differences that arise from different formulations of disordered models on a lattice. For 
concreteness, let us consider another disordered model: 

Z [{f3}, n] = Tr Si J] exp {(3 ijSi ■ Sj ) , (8) 

(ij) 

with a bond distribution P(f3). The formal operations lead to, again, a coupled model: 

M / M M 

(z[{(3},n]) M -U^ II ex p aE^-^+a £ ^M a M b M' 

a=l (ij) \ a=l a =l,b=l 

+ (higher order terms) , 

where (3 X = and fh = \ — A/) are the first and the second cumulants obtained 
from the distribution P(/3^), respectively. Comparing this model (9) with the previous 
one (7), we notice that the differences between these two are simply the presence of the 
exponential function and that of the diagonal coupling: /3 2 (sf 1 ^ ■ s^ sf^ ■ s^^j with 

a = b in (9). Both features lead to the term (s\ a ^ ■ s^^) i.e. the process of the double 
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occupancy of the bonds on the same replica plane. Such a term can also be generated 
in the model (7) under block spin transformation. However, when the corresponding 
pure model is in the dilute phase, the process is irrelevant. For this reason, we expect 
both of the coupled models (7) and (9) are described by a same field theory. 

A related model, which is more suitable for taking continuum limit, can be 
formulated with a probability distribution P s {fx) for disorders fx on each site: 

Z [{//}, n] = Tr Si Yl ex P i^i + A*i) s i ' s j ■ ( 10 ) 

(ij) 

In this model, the bond strength j3ij = fii + fij respects a distribution P obtained by 
the convolution of two P s i.e. P s * P s = P. If this distribution P(f3) is identical to 
the distribution P((3) in (8), then this model becomes similar to that of (8). Still, the 
model (10) is different from the model (8), since there exists the correlation between 
the nearest-neighbor bonds (say and f3 ik ) connected at a site. Nevertheless, the 
correlation is still short range; we expect a similar behavior at large scales. 

We can argue that this correlation existing in the model (10) is not essential to 
distinguish it from the model (8) by examining how these two models are mapped into 
the theory space of the coupled models. Averaging over the disorder fi in (10) yields, at a 
site i, a term J2 a bjk ' s< f^ s f^ ' as a leading nonlinear coupling. In particular, 
the j ^ k part of this coupling is induced by the correlation between the nearest-neighbor 
bonds in the model (10). If this coupling were, with some special reason, never generated 
in the model (9) mapped from (8), we would seriously distinguish (10) from (8) in the 
first place. But, it has already been contained in (9), implicitly. Indeed, although the 
expression (9) is written with respect to links, this can be recasted in terms of sites 
thanks to the discrete rotational and the discrete translational invariance of the coupled 
model on a lattice. Then we recognize, by expanding the exponential, this coupling is 
also present in the model (9). 

To summarize, in the dilute phase, there are two important processes on a 
lattice: the walk on the same layer (^s^ ■ s^^j and the branching to another layer 

(s^ ' s ^ s f^ ' f° r a b. It is plausible to think that these coupled models 
discussed here are, at large scales, described by a simple continuum field theory. 

2.2. Critical point of the homogeneous 0{n) model 

When we approach a critical point of some lattice model defined in the dimension d, 
the correlation length £ of the model, which is usually the order of the lattice cut-off 
a, becomes arbitrary large. The continuum limit (or, scaling limit) of the lattice model 
can be reached by taking a — > 0, while keeping £ fixed. There is nice class of local lattice 
observables {0^ at } which have the finite limit 

lima -2(A, 1+ ... + A,J (0 lat (;ri) ... lat (Xn)) = (^(^ . . . 0^)), (11) 
a—>0 

for an appropriate choice of scaling dimensions {2A^ fc }. With the idea of the block spin 
transformation, the scaling fields <f>k(xk) ma Y be considered as coarse grained versions 
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Lattice 



Continuum primary fields scaling dimension 



spin 
polarization 
encounter 



energy 




S(x) 0i,3 -2 + 4a 2 

a(x) p _ ljP 1 - a 2 _/2 - 3/(8a 2 ) 

a u (x) 0i jO 1 - ai /2 

<Jw{x) 2 ,o 1 - a 2 _/2 + 3/(2a 2 ) 



Table 1. Correspondence between lattice observables and continuum fields. The value 
of p is given by p = 1/ (2 — 2a 2 _ ) . 



of fc at over some region of size L, centered at Xk, with a « I « (. These scaling 
dimensions are determined dynamically from the interactions in the system, and related 
to the RG eigenvalues y^ via y^ = d — 2A^ fc . When y$ h is positive (negative), the RG 
flow near the critical point is unstable (stable), and the corresponding observable and 
fields are relevant (irrelevant). 

Around the critical point in the theory space of the 0{n) model, there are, at least, 
two directions which are unstable {i.e. relevant) under the RG flow to the infrared. 
These directions are spanned by the coupling constants associated with the two most 
important pairs {0|f t , 0aJ of local lattice observables and the corresponding scaling fields. 
On the lattice side, one is the energy density s$ • Sj and the other is the spin vector 
sf. In the scaling limit, they renormalize to become the scaling fields called the energy 
operator £{x) and the spin operator a(x). The scaling dimension of £{r) and that of 
a(r) is related to the leading thermal and the leading magnetic eigenvalue of the spin 
system, respectively. These two positive eigenvalues correspond to the linearized RG 
flow near the critical point, and determine the thermodynamic exponents [16]. 

So far, we have described the qualitative structure of the theory space focusing on 
the correspondence between the lattice observables and the continuum scaling fields. It 
is tantalizing that there is no versatile scheme for a given lattice model to establish this 
type of correspondence and to proceed to quantitative analysis using a scaling theory at 
hand. However, the 0(n) model was particularly successful in this respect; the geometric 
nature such as the loop representation of the partition function made it possible to access 
some of the exact results on the scaling dimensions by the CG method [18] before the 
emergence of the conformal field theory [3]. These exact results by the CG method, or 
by the other means [14] then led to the conjecture claiming that a certain one-parameter 
family of CFT's describe the critical points of 0(n) model for the continuous value of 
n (\n\ < 2) [12]. This claim is further checked by the Bethe ansatz result [2]. 

We assume this one-parameter family of CFT's [12] in the rest of the paper. The 
correlation functions of primary fields are represented in terms of the vertex operators 
^ e ia k <t>{x) j an appropriate choice of charges {a k }, where <f>(x) is a bosonic scalar field 
[12, 15, 33]. The necessary formulae on the vertex operator representation are given in 
Appendix A. In this construction, we have two marginal operators (the screening vertex 
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operators). In the 0{n) model, the screening charge a_ is related to n as 

n = -2 cos (n/a 2 _) . (12) 

The energy operator and the spin operator is identified with the primary fields as follows: 

£(x) -> 0i i3 , cr(x) -> (pp-l, p , (13) 

where the value of p is given by p = 1/(2 — 2a 2 ). These identifications are based on the 
dimensions originally obtained by the CG methods given in Table 1. Using the formulas 
(A. 5) and (A. 6), the dimensions of these operators are checked as 

2A £ = -2a lj3 a I ^ = -2(-c^_)(a+ + 2a_) = -2 + 4a 2 , (14) 
2A CT = -2a p _i iP «^ T ^ = -- ((2 - p)a + + (1 - p)a_) (pa+ + (1 + p)ai_) 

= l-a 2 /2-3/(8a 2 ). (15) 

where = — 1 in (A. 3) is used. 

The relation to the two dimensional critical statistical models are summarized 
below. The 0(n) model at n = 2 (a 2 = 1) and n — 1 (a 2 . = §) are the XY model 
and the Ising model, respectively. The Ising model is the only minimal unitary model 
[3] that belongs to the 0(n) family. Note that the Ising model is also regarded as the 
g-state Potts model at q — 2 [8]. The model at n = (a 2 = |) corresponds to the 
polymer, or self avoiding walk (SAW), as the qualitative result was obtained by e = 4 — d 
expansion [26]. The model in two dimensions is related, under the SLE duality [32], to 
the percolation (Potts model at q — 1). Both models are non-unitary, and various 
scaling dimensions in the models are numerically studied in ref . [28] . Finally, the model 
at n = — 2 (a 2 = |) is the loop-erased random walk (LERW). The one-parameter family 
of 0(n) CFT's are related to the SLE by 

/t = 4a 2 , (16) 

where k is the strength of the Brownian motion which drives the evolution SLE. These 
critical 0{n) models covers 2 < k < 4. 

In the rest of the sub-section, we discuss more on the correspondence between 
the lattice and the continuum. The relation between the loop configuration and the 
correlation function in the continuum theory is given; the latter is the object of our 
study henceforth. The reader who is not interested in the lattice may skip the following. 

In this paper, we restrict our consideration to the calculation of the RG flow in 
the subspace spanned by the coupling constant associated with the energy operator £ 
and the spin operator a. There is, however, another important RG eigenvalue apart 
from the ones associated with these two; the eigenvalue is responsible for the geometric 
property of the loops. In other words, the thermal and the magnetic eigenvalues are 
not sufficient to characterize the local shape of the loops even at a primitive level. To 
see this, let us first remind that, in the SAW (n = 0), the fractal dimension D F of the 
loop is given by the thermal eigenvalue ys, and thus D F = \jv [26]; it is then given 
by ye = 2 — 2A^ = | in two dimensions. The generalization of this result to arbitrary 
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n {\n\ < 2) was discussed by relating the 0{n) loops to the clusters in the percolation 
[57]. The fractal dimension for n ^ should then be written as 

D F = —, (17) 

with 0> ^ 1 (in the percolation theory, a> is called the Fisher exponent, and v 
the correlation length exponent). Further, the fractal dimension Dp is expected 
to be the RG eigenvalue from the polarization operator (or, the two-leg operator): 
{0 lat , 0} = {sf sf , crn(x)} with the vector indexes a ^ [3, for which the scaling dimension 
can be obtained by the CG method (see Table 1). The fractal dimension of the 0(n) 
loop in the dilute phase is then given by D F = 2 — 2A CTn = 1 + a 2 _/2 + . 

Diagrammatically, the polarization operator an(x), when inserted into a correlation 
function, creates two curves emanating from a point x, while the spin operator a(x) 
creates only one curve. The curves created by o~u{x) repel each other; they have different 
colors a ^ (3 and hence can not connect by themselves. They need to be annihilated by 
the other operator au(y). By contrast, the insertion of the energy operator S(x) require 
that a point x be passed by some loop; the curve passing the point will connect each 
end to form the loop. 

However, forming a loop becomes more and more difficult as n tends to zero; the 
loop segments repel each other in the SAW. In this sense, the role of the energy operator 
8 approaches that of the polarization operator an. In fact, at the SAW, these two 
operators has the same scaling dimension 2A = |, and this explain ap — 1 in (17). In 
the CFT context, by using the scaling dimensions obtained by the CG method, an(x) 
was argued to correspond to the 01,0 primary field. Then, ap = 1 can be confirmed by 
the well-known equivalence between the primary field: 

4>r,s = 4>m-r, m +l-s (0 < T < 171, < S < Ul + 1) , (18) 

in the m-th minimal model with the central charge c = 1 — 6/m(m + 1). Indeed, 
in the SAW, which is the m = 2 minimal model, 0^3 (the energy operator £) 
should have the same dimension as the field 01,0 (the polarization operator an), since 
(1,3) + (1,0) = (2,3). 

There are, of course, many other composite fields constructed as products of the 
several local lattice observables. A quantitative discussion on such fields is given, which 
is based on the OPE and the correlation inequalities [19]. Their result suggests, in 
general situation, that the scaling dimension of a composite field exceeds the sum of 
those of the elementary fields due to the repulsion between them. This implies higher- 
order composite fields are more likely to be irrelevant in the RG. In our context, we 
mention the four-leg operator {J2j( s i ' s j) 2 > ^viO^)} as an important example [18, 33]. 
The insertion of the field ayi(x) into a correlation function corresponds to requiring two 
of the loop segments to overlap; in the particle picture, the trajectories should encounter 
at the point x. Again, the dimension of the field ayi(x) is obtained by the CG method, 
and identified with the primary field 02,o in the CFT. Now, we see that it is relevant in 

+ This fractal dimension is also derived, in a rigorous way, from the SLE [32] 
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the dense phase (k > 4) and irrelevant in the dilute phase (k < 4) as shown in Fig 2. 
As used in Section 2.1, the latter observation is a key ingredient in understanding the 
universality of the pure 0(n) model in the dilute phase. Further, at the level of coupled 
model, this is crucial in our expectation that (7) and (9) are, near the decoupled critical 
points, described by a same field theory. 

2.3. Scaling limit of the disordered 0(n) models 

We shall go into the continuum formulation of the disordered models by using the CFT 
description of the critical 0(n) model. In section 2.1, we have considered, on the lattice, 
three disordered models whose partition functions are given by (3), (8) and (10). As we 
have seen, after the disorder averaging, they have much in common when mapped to the 
coupled models. When the couplings between replicas are sufficiently small, one could, 
by looking at (7) for instance, guess the action of the coupled models in the continuum 
limit. But we shall take the other way. Instead, we first consider the continuum limit of 
the disordered model (10), and then take the disorder average in the continuum theory. 
The advantage is that we can use the OPE to deal with the composite operators. 
In the following, we will use the relation: 



1 +tsi ■ Sj) ~ exp 



■,3) 



exp 



(19) 



The first equivalence symbol means that these two pure model belong to the same 
universality class, as discussed in the paragraph above (8). In the right hand side, 
we formally write the action of the 0(n) CFT as S*. The coupling constant m is 
proportional to the reduced temperature, which is given by (t — t c )/t c , or (f3 — f3 c )//3 c . 

Using the relation (19), it is natural to assume the continuum limit of the model 
(10) is described by 




,(20) 



(Z[{^},n]) M = I f[V& a \x) exp - ) + f d 2 xm(x) 

J a=l |_ \a=l / •* 

where Si a ^ is the action of the 0(n) CFT on the replica plane (a). Here, we formally 
write elementary fields (in the sense of [7]) as & a \x) and replace the tracing Tr <„) in 

(10) to J £>$ a (x), which denotes the path integration over the fluctuation of & a \x). 
The realization {//} of the local weights //j on lattice sites in (10) are now replaced to its 
continuum counterpart, that is, a configuration of a scalar function m(x). Physically, 
this m(x) serves as a locally fluctuating reduced temperature. We assume this m(x) 
respects a single distribution function P (m(x)) which is independent of the position 
x. Hereafter, we suppress the argument of the partition function Z [{//}, n]. Then the 
averaging operation is written as 

J Vm(x)P{m(x)) Z M , (21) 



Z M 



Disordered 0(n) Loop Model and Coupled Conformal Field Theories 



15 




Figure 2. The scaling dimensions 2A of the primary fields ^3, <f>\^ and ^2,0 plotted 
as functions of the SLE parameter k = 4a_. The horizontal line at 2 A = 1 shows the 
marginality of the most relevant couping £^(x)£^(x). 



where T>m(x) denotes a path integral measure Y\_ x dm{x). Let be a fc-th cumulant 
of the distribution function P(m(x)). In the important case of a Gaussian distribution 
P(m(x)) = exp [—(m(x) — m ) 2 /2<?o]; we have £1 = m , £ 2 = 9o an d Ck — for k > 3. 
Then, by the averaging over the disorder, we obtain 



M 



(22) 



Z" = f f[V<5> a (x) exp - J2 S i a) + [ d2x Ui^S^x) + Sj( 

J a=l I a=l V o=l 

M M 

Si(x) = i 2 Y,£ {a \x)£ ib \x) + ^Y. £[a) ^ £ib) ^ £iC) ^ + ■■■■ ( 23 ) 

a, b a,b,c 

In this formal expression, we should note that the nonlinear part Si(x) contains the 
composite operators i.e. the products of the several fields at the same point on the 
same replica plane. We deal with them by using the fusion rules in the 0(n) model. 
According to the identification £ — > 13 in (13), the fusion rule in the thermal sector 
reads 

01,3 • 01,3 ~ 01,1 + 01,3 + 01,5, (24) 

or, equivalent ly, 

£ ■ £ ~ / + £ + £'. (25) 

Here, £' = 0is is the next-leading energy operator. Since this operator is irrelevant 
as shown in Figure 2, we neglect the £ ■ £ — > £' part of the OPE (25). The identity 
field I = 0n contributes as a trivial shift of the free energy. The emergence of the 
energy operator £ in the right hand side of (25) is a characteristic of 0(n) model 
with n 7^ 1, which will be discussed in the paragraph below (35). Because of this 
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contribution (£•£—>£), qualitatively speaking, we should have a hierarchical flow of 
the coupling constants: • • • — > £3 — > £2 — ¥ ii- For example, the diagonal (a = b) part of 
£2 Yl £ i - a \x)£ <yh \x) mixes with £1 J2£^ and hence the flow £ 2 — ► £1 occurs *. Thus, by 
redefining the coupling constants , we get 



M 



m / M 



(26) 



Z M = ff[V<Z> a (x) exp - J2 S * a)+ [ d2x [m J2 S{a) ( x ) + S i( x ) 

J a=l I a=l \ a=l 

M M 

Sj(x) = g Y,£ {a) (x)£ (b) (x) +6 £ {a) (x)£ (b) (x)£ (c) (x) + ■■■, (27) 

a^b a^b,bj^c,Cy^a 

where we use the symbol m Q , g and ^ (k > 3) for the new coupling constants. Now, the 
effective interaction Si(x) is defined without composite operators. It should be noted 
that we restrict our considerations on the theory space which is replica symmetric; we 
assume that the coupling constant for, say £^ (x)£ ( - b \x), is independent of the pair 
(a, b). 

These coupling constants are determined by the distribution function P(m(x)) in 
(21), and are non-zero in general. The massless limit (m — > 0) of the decoupled model 
(g = 0, £fc = 0) remains obviously as a fixed point. We now change the distribution 
function P(m(x)), and gradually turn on the effective interaction SV in (27) while 
keeping m = 0. The large scale behavior is then dominated by the most relevant 
field: £^(x)£ ( - b \x). The dimension of this operator is given by 4Ag, which is the twice 
the dimension of the energy operator. 

As we can infer from (14), this field becomes marginal at n — 1 (a 2 = |). Since 
0(1) model is the Ising model, this serves as a simple check of the known marginality of 
the disorder in the random-bond Ising model [8, 30]. Hence, we use the parametrization: 

a - = \~ e ' ( 28 ) 
to perform the epsilon expansion in the next section. 



3. Renormalization group flow in the one-loop calculation 

In this section, we discuss the scale dependence of the theory (27) by one-loop RG 
calculation. It will turn out that under certain circumstances, the disordered 0{n) 
model has a pair of the ultraviolet (UV) and the infrared (IR) fixed points, as in the 
random-bond g-state Potts model [8, 20]. We will investigate the properties of the IR 
fixed point by which the large scale behavior of the theory is dominated. 



3.1. The epsilon expansion and procedure of RG 

We shall calculate the RG flow of the coupled model (27) perturbatively in the epsilon 
expansion using the parameter in (28). Generally, a flow is determined only by the 

* In the coupled model at M ^ 0, this flow causes shift in the critical temperature. The magnitude of 
this effect is proportional to the number of diagonal elements M, the structure constant (35) and £2- 
However, since the disordered model corresponds to M — > 0, we assume it is negligible. 



Disordered 0{n) Loop Model and Coupled Conformal Field Theories 



17 



most fundamental properties of a theory such as the symmetries and the dimension of 
the space. In order to grasp the topology of a theory space, it is often useful to think 
that these symmetries are dependent on some continuous parameters [7, 8, 20, 24], or 
that the dimension itself is a continuous parameter [25] . Although the change is gradual 
when we look locally at a generic point as varying these parameters, the global structure 
of the flow, on the other hand, may change drastically at certain critical values of the 
parameters. A typical topological transition is caused by a branching of one merged 
fixed point into a pair of the UV fixed point and the IR fixed point. The idea of epsilon 
expansion provides us with a firm ground to perform a perturbation calculation in the 
non-trivial region after the branching. 

In the parameter space, we can infer the location of the branching by the presence 
of a marginal field: the hallmark of the merged fixed points. In the renowned example 
of the ^-invariant scalar field theory [25], the most relevant interaction 4 becomes 
marginal at the dimension d — 4. In d = 4 — e, the Gaussian fixed point and the 
Wilson-Fisher fixed point emerge; we can do a solid perturbative calculation by setting 
an appropriate coupling coordinate in which the coupling constant remains 0(e). 

In the disordered 0{n) model, as mentioned in the previous section, the most 
relevant part of the effective interaction Sj in (27) is the term £ < - a \x)£( b \x), which is 
marginal at n — 1 and is relevant for n < 1 (see Fig. 2). The next-leading relevant 
term £ (a \x)£ {b \x)£^ (x) is irrelevant for n > and is marginal at n = (SAW) J. We 
restrict our analysis on n > 0, keep the most relevant part 



and discard the other higher-order terms in (27). 

The implementation of the RG is as follows. We introduce a cut-off length scale r, 
which serves as an effective lattice spacing, and determine a renormalized coupling g(r) 
by integrating out the short-distance degrees of freedom: 



r- 8 ^(r)(^(x)^(oo)) =g {S I (x)S I (oo)) + ^ / €py{S I (x)S I (y)S I (oo)) 



where the insertions of the interaction fields Si are restricted onto a disk of radius r, 
and the coupling is measured by projecting perturbative contributions on Sj(oo). The 
symbol (...)o represents the unperturbed correlation function. A trivial scaling factor 
r~ 8e is introduced in order to make the renormalized coupling dimensionless. At this 
stage in RG, as is well known, a finite redefinition of the coupling is possible; physical 
quantities (scaling dimensions, for instance) is invariant under a finite coordinate 

| It may be noted that the role of this field is analogous to that of the <j> 6 term in the Z2-invariant 
scalar field theory. This term is irrelevant if d > 3 and becomes marginal at d = 3. In 3 < d < 4, we 
have only the two fixed points mentioned above. 



M 




(29) 




(30) 
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Figure 3. The diagrams for the beta function. A gray disk and a nearby pair of the 
disks represent an energy operator J2 a £^ and the interaction Si = X)a/b^ a ^^' 
respectively. An arrow corresponds to the sub-leading part of the OPE in (34) and a 
shaded-square represents four-point function. The external lines to the point infinity 
in (30) are amputated. The second order: (a) G2,i(r,e) and (b) G2,2(r, e); the third 
order: (c) G 3 ,i(r,e), (d) G 3 , 2 (r,e), (e) G 3;3 (r,e) and (f) G 3 , 4 (r,e). 



reparametrization of theory space [24]. To fix this ambiguity, we choose a minimal 
subtraction scheme. Then we get the renormalized coupling in the form: 

g(r) = r 8£ (g + G 2 (e, r)g 2 + G 3 (e, r)g 3 + •••), (31) 

where Gz(e, r) and Gs(e, r) have only poles in e and no regular part §. The beta function 
is then obtained by differentiating the renormalized coupling (31) with respect to (lnr). 
Using the first two terms in the fusion rule (25), we make contractions and list, in Figure 
3, the possible diagrams for the beta function. 

3.2. One-loop beta function, a line of random fixed points and a strong coupling region 

We here calculate one-loop beta function. The process represented by the diagram in 
Figure 3- (a) involves three layers, and the number of the ways for this contraction is 
4(M - 2). Then we have 



G 2 ,i(r, e) = ± • 4(M - 2) / d 2 y(£(x)£(y)), 

z - J\v-x\<r 



'\y-x\<r 

47r (M-2)(^V (32) 



-4Af 



where we have used (£(x)£(y)) — \x — y\~ £ and 4A^- = 2 — 8e. 

Besides this, up to one-loop order, there is the other contribution represented by 
Figure 3-(b) which is obtained by using twice the sub-leading part of the OPE {£■£ — > £): 



G 2i2 (r,e) = -.2 / dy f ] — — — 1 

l - J\y-x\<r \\ X V\ J 



y—x\<r 

§ More precisely, we keep the combination r 8£ /8e = l/8e + logr + (9(e) and discard the other parts. 
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= 2*[C§ £ (c£)] 2 (^y (33) 

Here, C ££ (a 2 _) is the structure constant of the CFT which appear in the OPE: 

£ (0, 0) . £(,, *) = A, + _i&_ £ + ^ + . . , (34) 

where we write just the primary operators to represent their conformal families including 
the descendants. The structure constants are, in general, determined by the requirement 
that the operator algebra be associative. In practice, the crossing symmetry of the four- 
point function is strong enough to fix them [3]; the actual values are obtained by using 
either the connection matrix of the hypergeometric function [11], or the vertex operator 
representation of the four-point function [13]. This reads, 

C e M = 2(1 - 2p) ^ 7(3 _ 4p) , (35) 

where we have used temporary notation p = a 2 _ and 7(2;) = T(x)/T{\ — x). 

The structure constant C ££ {a 2 _) has information about the important selection rules 
in the pure 0(n) model. First, observe that the square of the structure constant (35) 
has a zero at a 2 _ — | (the Ising model, e = 0) and a pole at a 2 _ = | (SAW, e = ^). 
The former zero is due to the self-duality of the Ising model [17]. In the vicinity of the 
critical point, the duality transformation changes the sign of the reduced temperature 
t. Since the energy operator £ couples to t, if the model is invariant under the duality, 
it should be odd under the duality. Then £ is not allowed to appear in the right hand 
side of the fusion rule (25). 

The pole at a 2 _ — | emerges because of the strong repulsion between the loop 
segments in the limit n — > 0. Since the normalization of operators are fixed such that 
C ££ = 1, what really happens is the divergence of the ratio C ££ /C ££ . When the two 
loop segments approach, they will either (i) form a complete loop (to contribute the 
free energy) or (ii) form together a joined loop segment. Since, in the limit n — ■> the 
process (i) is strongly suppressed by the repulsion between the loop segments, the ratio, 
indeed, diverges. 

Taking these interpretations into the account, more intuitive representations of 
the second-order contributions are possible. The term G2,i(r,e) is represented by the 
diagram in which two open segments lie on two layers and one closed loop on another 
layer (Figure 4-(a)), while the term G 2: 2( r ,^) is represented by the diagram in which 
two open segments lie on two layers (Figure 4-(b)). 

Now we sum up (32) and (33) to get 

9(r) = r 8£ {so + 2vr [2(M - 2) + [C £ ee {a 2 _)] 2 ] (j^J g 2 + O(^ 3 )} . (36) 
By introducing g = r 8e g , we solve (36) with respect to g: 

g = g(r) - — [ 2 (M - 2) + [C £ ee {a 2 _)] 2 \ g(r) 2 + O (g{rf) . (37) 
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(a) 



(b) 







Figure 4. The second order diagrams; see also Figure 3. (a) G2,i(r, e): the three-layer 
process with a closed loop (b) G 2! 2( r > e ) : the two-layer process. 



Using this, the beta function up to one-loop order is then obtained as 




(38) 



This beta function of the coupled model permits an IR fixed point when the coefficient 
of g(r) 2 is negative. As we take the limit M — > for the disordered model, there 
are a competition between the terms from G 2) i(r,e) and G2,2(r, e). As a result, we get 
two distinct regions for < n < 1: in the first case we have a monotonic increasing 
beta function; the theory goes to the strong coupling. In the second case we have an 
IR fixed point. The one-loop beta function suggests that the IR fixed point vanishes 
at e = • 0.770861 • • •, which corresponds to n = n* = 0.26168 ■ • •. The flow of the 
coupling constant is schematically shown in Fig. 5. It should be noted, however, that the 
divergence of the coupling g at n = is an artifact of the one-loop calculation; if there 
is a positive g 3 term^f, the line of the IR fixed point terminates at (n c ,g c ) = (n* ,g' CJf ) 
with a finite coupling g' CSf . 

The strong coupling phase of the disordered model (M — > 0) near n m is formed 
by the dominance of the two-layer process (^2,2 ( r , e) over the three-layer process with a 
closed loop G ? 2,i(^, e)- The qualitative behavior of the absolute ratio lZ(n) = IG^/G^i] 
is independent of M except a special case with M = 2, where the three- layer process 
is prohibited. Diagrammatically, a large TZ(n) implies that the ratio of (the number of 
inter-layer hopping) to (the average number of complete loops per layer) is also large. 
Near n ^ region, the particles, which are destined to form a whole connected diagram, 
favor to sew layers together rather than to stay on one layer and to walk around avoiding 
their traces. 

We shall henceforth proceed to the two loop calculation to investigate the IR fixed 
points in the region near n = 1 (e ^), where regarding the structure constant 
in (35) as 0(e) is justifiable. In this region, the calculation turns out to be parallel to 
that in the study of the random-bond Potts model by Dotsenko, Picco and Pujol [20]. 
At one-loop level, Cf g = 0(e) implies (^2,2 = 0(e); in the minimal subtraction scheme, 

|| A similar form of beta function appears in e = d— 2 expansion of the random-bond Ising model [16]. 

Although the e dependence of the third order coefficient is not calculated in this paper, we anticipate 
a positive g 3 term from a continuation of the result at e <C in Section 4. 
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Figure 5. A schematic flow diagram of the couping constant g with respect to n of 
the disordered 0(n) model suggested by the one-loop beta function (38). The bold 
curve represents the positions of the IR fixed points. 



we drop this term. The beta function is then 

P{g) = = Seg + 4vr(M - 2)g 2 + 0(g 3 ). (39) 

We left the two loop calculation of the intermediate region where C^ £ « 1 as a future 
problem. 

4. The scaling dimensions up to the two-loop calculation 

4-1. The beta function at two-loop 

There are four types of third order diagrams for the beta function as listed in Figure 
3-(c), (d), (e) and (f). As the diagram G 3 , 3 (r, e) in Figure 3-(e) have extra e 2 factor from 
the structure constants and the diagram G 3i4 (r, e) in Figure 3-(f) has no pole in e, we 
omit these terms in the following. The G 3i i(r, e) in Figure 3-(c) can be calculated as 



G 3l i(r,e) = i-12(M-2)(M-3) / d 2 yd 2 z(£(x)£(y)) (£(y)£(z)) 

J \y—x\<r,\z~x\<r 

= 8tt(M-2)(M-3) / dyy~ 1+We [ d 

J\v\<r ^l<H 



' 2 z|d- 2+8e U-i|- 2+8e 



r 16e 

= 16vr 2 (M - 2)(M - 3)^ + 0(e°), (40) 
where we have used the asymptotic behavior of the integral for 1 <C R: 

d 2 z\z\- 2+ ^\z - l|- 2+8e = 4tt1 + 0(R- 2+16e ). (41) 

l\z\<R ° e 

It should be emphasized that we get no simple pole in (40). 

Next, we have G 3 ,2(?", e) in Figure 3-(d) containing four-point functions: 

G 3i2 (r, c ) = I • 24(M -2) I d 2 yd 2 z (£(x)£(y)£(z)£(oc)) Q (£(y)£(z)) Q (42) 

^/ |y — a;|<r,|^; — a;|<r 



/ 
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J y<r 

where we define 



l(R,e) 



\<R 



d 2 z(S(0)S(l)S(z)S(oo)) {£(l)£(z)) . 



22 
(43) 

(44) 



Actually, instead of X(R,e), we shall calculate Z(oo,e) by analytic continuation in 
Appendix B, and see that Z(oo, e) itself has no poles in e. Although, in the Appendix, we 
prove this fact by applying the contiguity relation between generalized hypergeometric 
series, the fact can also be interpreted as a result of a cancellation between two poles 
with opposite sign, as explained in the following + . 

In the limit \z\ — > oo, the integrand in (44) behaves as |z|~ 4A£ (4A £ = 2 — 8e), since 
the four-point function, using the identity operator / in (34) as the leading intermediate 
channel, decouples into two two-point functions. Hence, Z(oo, e) has the contributions 
from both regions \z\ < R and \z\ > R, each of which contains a pole +27ri? 8e /8e and 
— 27i\R 8e /8e, respectively; they are canceled out each other. Observing that I(R, e) does 
not have the latter pole, and estimating the corrections, we obtain 



8c 



T(R, e) = 2tt^- + J(oo, e) + 0(iT 1+8e ). 



(45) 



As the first term corresponds to the decoupling limit y — > x of the four-point function 
in (42), it leads to the same contribution as (40), only if we replace (M — 2) in (43) by 
(M — 2)(M — 3). Now, using the vertex operator representation, we write Z(oo,e) in 
(45) as 



J(oo,e)=AT J d 2 zd 2 ud 2 v (^3(0)^3(1)^(^3 



[oc)V a _(u)V a _(v))(V ai3 (l)V^(z)) 



d 2 zd 2 ud 2 v |d 4ai3ai3 |l - z \^ a ^\ 



'z — u)(z — v] 

1 4a 2 



\Aa13a- 



(46) 



• \uv(u- l)(v- l)\ 4ai3a - \u - v\ 

where M is the normalization factor with the four-point function given in (A. 7). 
Substituting (45) and the result (B.30) for Z(oo,e) into (43), we obtain 

1 1 



G 3 , 2 (r,e) = 16vr 2 (M-2) 



(47) 



r 16e + O(e ). 

Sty 8e_ 

We get, by collecting the third order terms (40) and (47), the renormalized coupling 
constant in (31) as 

'(M-2) 2 (M-2) 



g(r) 



„8c 



flo + 47r (M Q 2) r 8 ^ 2 + 167r 2 



8e 



r 16e 9 3 } + O(g*). (48) 



Then (48) is solved with respect to g — r 8e g as 



g = g (r) - 4tt 



(M - 2) 2 (M - 2) 
(8e) 2 + 8e 



g{rf 



(49) 



A similar cancellation of poles has already been observed in the random-bond Potts model [20] . 
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Consequently, we obtain the RG beta function 

P(g) = J^ = 8e 9 + 4tt(M - 2)g 2 - 16vr 2 (M - 2)^ 3 + 0(g 4 ). (50) 

As is expected from the known equivalence of the random-bond Ising model to the 
random- mass fermion model and the Gross- Neveu model [30, 8], at the Ising point e = 
(n = 1), the expression (50) reduces to the two-loop beta function for the Gross- Neveu 
model [34]. Solving P(g c ) = for the disordered model (M — > 0), we obtain the coupling 
constant at the random fixed point as 

g c = - + — + 0(e 3 ). (51) 

71 71 

4-2. The correction to the scaling dimensions 

In our perturbation theory, the two-point correlation function between fields O(0) and 
O(oo) is calculated as 

(O(0)O(oo)> = (0(O)0(oo)) o + go I d 2 x (^(y)C(0)C»(oo))o 

JT ■ (52) 



+ % [ d 2 yd 2 z (^(y)^(^)O(0)O(oo)) + • • • 
J\y\< r , \ z \< r 



Again, as in (30), we restrict the insertion of the interaction Si around the field 0(0) 
within the disk of radius r. From this, we can determine a renormalization constant Zq 
for the field O perturbatively as 

(O(0)O(oo)) = Z o (O(0)O(oo)) . (53) 

The two-point functions in the theories with different cut-offs r and sr are related as 

(O(0)O(sR)) srMsr) = ^^y ^ A °(O(0)O(R)) rMr} 



cxp 



2 /'" (Sr) ^7^ s-' A o(cm)0(ff)) r (5i) 

g(r) P{9) 



where, in the second line, we have introduced the anomalous dimension 

dZ (g) f . 

lM = -dfrT- (55) 

Now consider the large s behavior of (54) and let the upper limit g(rs) of the integral 
tend to g c of the random fixed point; the beta function (3(g) tends to zero, while the 
anomalous dimensions 7o((?), as we will see, remain finite both for the energy operator 
£ and the spin operator a. This means the integral is dominated by the contribution 
from the region g rs g c , and we obtain 

2(A I R -Ag v ) =- 7 o(ffc), (56) 

where we mean, by 2A^ and 2Ag v , the scaling dimension of the field O at the IR (the 
random) fixed point and the UV fixed point, respectively. 
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Figure 6. The contributions for the scaling dimensions. The diagrams (a), (b) (c) and 
(d) are for Zg, while (e) 1S2 (r, e) and (f) Ss(r, e) are for Z a . The white circle represents 
the spin operator a^; for symbols, see also the caption in Figure 3. 



4-3. Scaling dimension of the energy operator 



The first and second order diagrams for the renormalization constant of the energy 
operator Z £ are listed in Figure 6- (a), (b), (c) and (d). As the integral in the diagram in 
Figure 6-(c) does not have pole in e, we omit this term from the calculation. Since the 
integrals are same as those for the two loop beta function, by adapting combinatorial 
factors, we have 



Z £ = 1 + g 



2(M 



2(M 



(M 
= 1 + 4wg v 



8e 



^2)° 2 > l + 2! 
^ + 8n 2 g 2 



\(M-1)(M 



2) 4(M 

^3,1 + 



2(M-2)(M-3)"°' 1 ' 4(M-2) 
(M-l)(2M-3) (M-l 



G 



3,2 



(57) 

5er 16e 

where the results (32), (40) and (47) are used in the second line. Again, by using the 
expression for g in (49), we get 

d In Zs 



ie{g) 



d lnr 



4tt(M - l)g(r) - 8tt 2 (M - l)g{r) 2 + O (g(r) 3 ) 



Consequently, by using (51) in (56), we obtain 

2 (Af - A^V) = - l£ ( gc ) = 4e + 8e 2 + G(e 3 ). 



(58) 



(59) 



4-4- Scaling dimension of the spin operator 

We shall here calculate the renormalization constant for spin operator Z a up to the third 
order in g , which gives the lowest order correction to the scaling dimension 2A CT at the 
IR fixed point. According to the operator algebra (see the diagrams in Figure 6), there 
are no contribution at the first order, and one contribution at the second order <S 2 (r, e) 
and the other one at the third order Ss(r, e): 



Z a = l+S 2 (r,e)+S 3 (r,e) + O(g 4 ) 



(60) 
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The second order diagram in Figure 6-(e) is given by 

2 r 

<S 2 (r,e) = |-4(M-l) / d 2 yd 2 z (a(x)£(y)£(z)a(oc)) (£(y)£(z)) , (61) 

J\y—x\, \z— x\<r 

while the third order diagram in Figure 6-(f) is given by 
Ss(r,e) = |f -24(M-l)(M-2) 

• / d 2 yd 2 zd 2 w(a(x)£(y)£(z)a(oo)) (£(y)£(w))o(£(w)£(z))o. (62) 

J \y—x\<r, \z— x\<r, \w— x\<r 

In evaluating the integrals (61) and (62), we keep only the terms with the lowest powers 
in e. For this reason, it turns out that we can add certain extra regions to these integrals. 
First, by adding the region \z — x\ > r to (61), we get 

S 2 (r,e)=Airg 2 (M-l) [ dy y- 1+1 ^ /C 2 (oo, e) , (63) 

J y<r 

with /C 2 (oo,e) defined as 

/C 2 (oo,e)= / d 2 z(a(0)£(l)£(z)a(oo)) (£(l)£(z)) . (64) 
Jc 

Second, if we add the region \w — x\ > r and \z — x\ > r to (62), we see that the third 
order contribution Ss(r, e) has the similar structure as the second order one <S 2 (r, e) in 
(61). In fact, under a trivial change of variables, we obtain the factorized form: 

<S 3 (r,e) =8ng 3 (M-l)(M-2) [ dy y~ l+2 ^ [ d 2 w\w(w - 1)|- 4A£ /C 3 (oo, e) , (65) 

Jy<r JC 

with the definition 

/C 3 (oo, e)= / d 2 z\l-z\ 2 ' 8A£ (a(0)£(l)£(z)a(oo)} . (66) 
Jc 

Then, we write the integrals (64) and (66) using the vertex operator representation: 
/C 2 (oo,e) =uj d 2 zd 2 ud 2 v |l-^|-^ £ (^a^(0)ya 1 3(l)K, 13 (^)^_ 1 , p (oo)y a _(«)y a _(7;)) 

■ \{l - u){\ - v){z - u){z - v)\ 4ai3a -, (67) 

/C 3 (00,e)=AT J d 2 zd 2 ud 2 V \ z f^=^^\\ _ z |4a 13 a 13 +2-8A £ |^|4a^a_| u _ v j4a 2 _ 

• |(1 - u)(l - u)(z -«)(*- v)| 4ai3a -. (68) 

and give the results in Appendix C. 

Now we obtain, by collecting the contributions (63) and (65), the renormalization 
constant for the spin operator as 



In Z a = Ang 2 Q {M - l)_/C 2 (oo, e) + 8ng 3 (M - 1)(M - 2) 



r 24, ^ 



24e 2e 



/C 3 (oo,e) + ^ 4 ), 
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where we have used (41) to obtain the third order term. Differentiating this with respect 
to (lnr) and using the result (C.10) for the integrals /C 2 (oo,e) and /C 3 (oo,e), we obtain 
the anomalous dimension for the spin operator as 

la {g) = 4vr(M - 1) [g(r) — (M — 2) <?(r) 2 } 2 • 2A/Ve [U 2 + W 2 - SYZne] 

+ 8tt(M - 1)(M - 2) (^j g{rf ■ 2A/Ve [U 2 + W 2 - 12FZ7re] + 0(g{r) 4 ) 

= 8A/V(M - 1) {eg{rf [U 2 + W 2 - 8YZ] - 4g(r) 3 (M - 2)YZ) + £>(#(r) 4 ),(69) 

where the expression for the bare coupling constant (49) is used in the second equality. 
For the definition of the constants A/", U, W, Y and Z, see (A. 7) and (C.5). The limit 
M — > for the disordered 0(n) model yields, 

la (g) = -8Afn 2 eg(r) 2 [U 2 + W 2 } + 0{g{rf). (70) 
Now, from (56), we obtain 

2(A IR_ A U V) = _ lM 



T 4 (-) 

'--Ye 3 + 0(e A ) 



where, in the last line, we have used the elliptic integral of the first kind defined by 

K(k) = f 2 i = ■ («) 
Jo V 1 - k 2 sin 2 

for the purpose of the comparison with the result in the random-bond Potts model [20]; 
this will be discussed in Section 6. 

5. Effective central charge 

The central charge, in general, characterize a system by counting the number of the 
critical fluctuations in it. If the replica method is used, however, the central charge 
becomes rather trivial; we are always left with the vanishing central charge c = 
because of the formal limit M — * in the end. Nevertheless, we can define a so-called 
effective central charge characterizing a random fixed point as follows [9] . Consider the 
M layers of the coupled model and the relation 

c IR (M) =Mc uv + Ac(M), (73) 

where cjjv an d cm(M) are the central charge of the pure system and that of the non- 
trivial fixed point, respectively. The effective central charge is then defined as 
dc m (M) 



'eff ~ 



dM 



(74) 

M=0 

In order to obtain the cir(M), it is useful to recall the definition of the C-function 
and the differential equation satisfied by it [22, 16]. To this end, we consider the vicinity 
of some critical theory S* in the theory space spanned by the coupling constants {g 1 }: 

S = S*- J d 2 x g%(x), (75) 
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and then specialize to our disordered 0(n) model. The response of the action under the 
transformation z^ — > z^ + a^(z) can be expressed by the stress tensor as 

5S = [ d 2 r T^a v . (76) 



2tt _ 

As usual, a particular component of the stress tensor T zz = |(T n — T 2 2 — 2iT 12 ) and 
the trace of it 4T ZZ = Tn + T22 are denoted as T and 0, respectively. As the RG flow 
dg l = f3 l dt occurs under the dilatation z^ — > (1 + dt)z^, the definition (75) and (76) 
leads to the relation: 

&(x) = 2n(3%(x). (77) 

From T and 6, one can consider three-types of the two-point functions: 

(T(z,z)T(0,0)) = F(t)/z\ (78) 

(T(z,z)e(0,0)) = H(r)/z 3 z, (79) 

(e(z,z)e(0,0)) = G(r)/z 2 z 2 , (80) 

which are measured at the scale r = ln(zz). The C-function is then defined as 

C = 2F-H-^G. (81) 

The conservation of the stress tensor d^T^ = yields d s T(z,z) + \d z Q(z ) z) = 0, and 
one can easily show 

i c -\ a < 82 » 

After fixing the reference scale r to r = 0, the C-function depends only on the coupling 
constants g l . We introduce the Zamolodchikov metric on the theory space as 

g ij = (zz) 2 ($ i (z,z)Q j (0,0)) . (83) 

zz=l 

From (77) and (82), we have 

\^^-C = -\{2,fg^^. (84) 

For unitary theories, the positivity condition implies that the metric should be positive 
definite. The C-function is stable at the point with (3 l = i.e. a fixed point and take 
the same value as the central charge as we can see from the definition (78)-(81). 
In our case, the perturbing field is quadratic in energy operator as in (29): 

$ = Sj = J2£ {a) ( x )£ (b) ( x ), (85) 

and hence the metric (83) is 

g = 2M(M -1). (86) 
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This metric becomes negative in the limit M — > 0, thereby violating the c-theorem. Note 
also that our normalization for the energy operator £ is such that C\ e = 1. Substituting 
the beta function (39) into (84), we obtain 

AC(M) = -6tt 2 • 2M(M - 1) ^ ? + C(e 4 ) 

1 ; 1 ; 6[4vr(M-2)] 2 V ; 

From the definition of the effective central charge (74), we have 

c S-^ V = 166 3 + 0(6 4 ). (88) 

Since we have the IR fixed points in e > 0, this always increase under the RG flow as 
expected. 



6. Discussion and Conclusions 



Before concluding this paper, we discuss the relation to the known results and then 
comment on possible future directions. 

The RG flow in our disordered 0(n) model near the disordered Ising point, which 
is inferred from the results (50), (58) and (69), is qualitatively similar to that in the 
random-bond g-state Potts model [8, 20]. The crossover to an IR fixed point occurs only 
when the most relevant interaction, which is quadratic in energy operators as in (85), 
becomes relevant; the region correspond to n < 1 for the disordered 0(n) model and 
q > 2 for the random-bond g-state Potts model. This can be understood by recalling 
how the two families of the pure models, namely, the 0(n) models and the g-state Potts 
models coincide at the Ising point (n — 1, q — 2) [12, 17]. 

The energy operator in the 0(n) model is identified with the primary field 0i i3 
as we have stated in (13), while the energy operator in the Potts model is identified 
with the field 2 ,i- These primary fields are, in a general minimal CFT, different 
objects as one expects from the fact that a four-point function of primary fields <p ryS 
is determined from an ordinary differential equation of order rs. Accordingly, in the 
vertex operator representation of the four-point functions, we should include two of the 
screening operators V- in the 0(n) model and one of the other screening operator V + 
in the Potts model, respectively. However, the field 0i j3 reduces to the other field 2 ,i 
at the Ising point (the m = 3 minimal model) because of the equivalence (18). At the 
level of the scaling dimensions, the reduction can be summarized in Figure 7, in which 
the line of 2A 1)3 and the curve of 2A 2 ,i intersect at the Ising point (k — 3). 

Up to one-loop order, when the deviation parameter e from the Ising point defined 
in (28) is small and neglecting the structure constant (35) in the 0(n) model can be 
justified, the beta function is determined only from the scaling dimension of the energy 
operator as in (39); thus similarity with the random-bond g-state Potts model is natural 
from Figure 7. By contrast, the reduction of the two-loop beta function (50), which 




Figure 7. The scaling dimension of the energy operators for the 0{n) model (2A 13 ) 
and for the g-state Potts model (2A 2 ,i) plotted as functions of the SLE parameter 
k = Aa 2 _. The horizontal line at 2Af = 1 shows the marginality of the disorder- 
induced coupling. The arrows indicate the directions of the expansion parameters e 
and ep. 



involves four-point functions, to that of the q = 2 random-bond Potts model in the 
limit n — > 1 serves as a non-trivial check of the vertex operator representations. 

Remarkably, two lines of the IR fixed points lie the opposite sides of the Ising point; 
if we introduce a parameter 

e K = 3 - «, (89) 

the random fixed points of the disordered 0(n) models and those of the random-bond 
g-state Potts models exist in the region e K > and e K < 0, respectively. Another natural 
parameter ep for the Potts models can be taken as a 2 + = |— ep, where the IR fixed points 
emerges in the region e P > *. The definition (28), the relation (16) and a 2 + a 2 _ = 1 
from (A. 3) lead to relations: 

1 9 27 

e= A eK = 'TQ ep ~6A e2p + 0{4) - (90) 
Then the results for the scaling dimension of the energy operator in the random fixed 
point of the disordered 0(n) model (59) and that of the random-bond g-state Potts 
models [8, 20] can be summarized in one expression as 

1 2 , „ f 1 + k + 0(el) (e K > 0; 0(n)) 

2Af = l + -(l-2A £ ) 2 + ((1 - 2A^) 3 ) = <^ 2 (91) 

[ 1 + -e 2 K + 0(el) (e K < 0; Potts), 

where both results are measured from the Ising value 2 As = 1. Naturally, the scaling 
dimensions at these IR fixed points imply faster decays in the correlation functions: 

* It should be noted that yet another parameter e defined as a 2 + = | + e is used in [20] . 
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Ag < A^ R as one expects with a general RG flow [27]. The ratio between the coefficients 
in e 2 . is |, which is the square of the ratio of the derivative of the line 2Ai ;3 at the Ising 
point to that of the curve 2A2,i in Figure 7. 

Contrastingly, the coefficients in the spin scaling dimensions are transcendental 
numbers. Using the relation (90), our result (71) for the disordered 0(n) model is 
written as 

2 (A? - A™) 0(b) = 2^J^k + O (4) (e K > 0). (92) 
Then we quote, from ref. [20], the result for the random-bond g-state Potts model: 

2 (A* R - A™)p ft = - 1 3j 9 U V p + 0(e 4 P ) (e P > 0) 
I Ws 32r(-I) 2 r(I) 2 (93) 

= d^[ K ( sin K ( cos Y2)] 2 ( - £k)3 + ° (e * < 0) ' 

where, in the second line, we have used the elliptic integral defined in (72). As a 
first observation, notice that both coefficients are consisting of the numbers K(cos9) 
known as "singular values" of the elliptic integrals for which the modular ratios 
K (cos 9) / K (sin 9) belong to quadratic irrational numbers. The ratios here are \fl 
and for the disordered 0(n) model and the random-bond g-state Potts model, 
respectively!!. These modular ratios are exceptional in that the modular angles 9 are 
commensurate to n, namely, 9 = tt/4 and 9 = 71/ 12. 

The coefficient in (92) and that in (93) characterize the universality classes of the 
disordered models as well as of the corresponding coupled two-dimensional CFT's in the 
zero-layer limits. Now, it is noteworthy that they are simply related to the quantities 
originating from the specific structures of three dimensional lattices; the coefficient in 
(92) for the disordered 0(n) model and that in (93) for the random-bond g-state Potts 
model correspond to the Watson's triple integrals [41, 42] which represent the special 
values of Green's functions on the body centered cubic (bcc) lattice and on the face 
centered cubic (fee) lattice, respectively. For our disordered 0(n) model, to be concrete, 
the coefficient in (92) can be written as 

2 K 2 (sinf) 1/1 r r r dudvdw 

7T 3 2n \ 7T 3 J Q J J 1 — COS U COS V COS w 

where the quantity in the parenthesis is the Watson integral for the bcc lattice. It 
would be interesting to see whether this coefficient can be derived, asymptotically at 
large scales, by a combinatoric analysis of coupled loop model on a lattice such as the 
one defined from the partition function (7). 

These remarks are also useful to give representations for the two coefficients in terms 
of rational integrals in order to see they are "periods" proposed in ref. [38] . Looking 
back the arguments, both coefficients are obtained not as products but as ratios between 

jj The two singular values have number theoretic origins and belong to the sequence of the values 
derived by Chowla-Selberg formula [40] ; the relevant part of the formula is obtained from the functional 
equations of the Epstein's zeta functions on the quadratic number fields Q(y/—T) and Q(y^3). 



(94) 
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the complex Selberg integrals [36] and the non-symmetric extensions of them such as the 
one obtained in Appendix D, and thus it is not a priory clear that they are periods even 
if each of these Selberg-type integrals independently turns out to be a period. But if 
we start from the Watson integrals, they can be transformed into rational integrals [41], 
and hence the coefficients in (92) and that in (93), characterizing the two universality 
classes interconnected at the disordered Ising CFT, belong to, in fact, periodsff. 
Now, we list future directions in the following: 

(i) The fractal dimensions and SLE with the replica symmetry. We have seen that 
there is a line of IR fixed points for n* < n < 1. An immediate issue to be discussed 
is the fractal dimension Dp of the loops in the IR fixed points. One may figure out 
the corrections for the Fisher exponent in (17) by perturbative calculation like the one 
performed in this paper. Related result has recently appeared for the random-bond q- 
state Potts model [62]. Further, the effective scale invariance on the line of the IR fixed 
points suggests the possibility of constructing one-parameter family of new SLE's. In 
order to describe the shape of the disordered loops, stochastic motion over the effective 
space enhanced by the replica permutation symmetry may be considered in analogy 
with the SLE driven by a composition of the stochastic motion in the real space and 
that in the internal space such as an SU(2) [50] or a Z N [51]. 

(ii) The strong coupling phase. We have seen that the coupled 0{n) model has 
moved to strongly-coupled theory by the two- layer process in Figure 4-(b) in the region 
< n < n*. The physical picture of the strong coupling region is important for the 
polymers in disordered environment [56]. Especially for this region, consideration on 
the theory space with the replica symmetry breaking (RSB) is desirable. Although the 
RSB-RG flow has been proposed in the random-bond g-state Potts model [21], numerical 
studies there supports the result of the replica symmetric fixed point rather than that 
of the RSB fixed point. We anticipate our model has more reason to favor the RSB 
situation due to formations of replica bound states boosted by the presence of the two- 
layer process. We should add that a large class of disordered models which includes 
strongly-coupled layered 0(n) models has recently been discussed from the view point 
of the AdS/CFT duality [61]. Actually, the focus of their analysis is on large-n, while 
our strong coupling region lies at n ~ 0; although a direct application seems to be 
difficult, this direction may be valuable in order to grasp the nature of strong coupling 
phases induced by disorder. 

(iii) Non-local observables. We have studied the scaling dimensions of the fields 
corresponding to the local segments of loops. In the pure 0{n) model, some results are 
known on non-local observables such as the area of loops [48] or the linking numbers 
around multiple points [49]. A challenging problem is the study of the off-diagonal 
pair-correlation between the global shapes of loops in the disordered model. The whole 
shape of a loop itself is, of course, difficult to handle. But extracting an essential 
part of such a correlation between non-local objects in a simple disordered system may 

ff More precisely, as suggested in the expression (94), the coefficient in (92) for the disordered 0(n) 
model belongs to the l/V-extended period ring V[1/tt] defined also in ref. [38]. 
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possibly provide us with the basis of other important problems. In chaotic systems, for 
instance, the off-diagonal correlation between periodic orbits, which are analogous to 
the disordered loops studied in this paper, is already a central issue in understanding 
the quantum level statistics [46, 47]. It would be interesting if the possible relation 
between the disordered loops and periodic orbits in a chaotic system could be clarified. 

In conclusion, we have studied the disordered 0(n) models in two dimensions. We 
formulated the inhomogeneous loop model on a lattice. After adopting the replica 
method in order to take the disorder average, the models were mapped to the coupled 
layers of the homogeneous loop models. Some of the possible differences in the lattice 
formulation of the disordered models were argued to be irrelevant using the properties 
of the pure model in the dilute phase. Then we discussed the continuum limit of the 
disordered model assuming the identification between the energy operator S and the 
primary field 01,3 in the vicinity of the pure critical 0{n) model. We considered the RG 
flow in the replica symmetric theory space. Since the most relevant interaction £^£^ 
becomes marginal at the disordered Ising model n — 1, we used the epsilon expansion 
methods near n — 1 to perform a perturbative calculation. At one loop level, the RG 
flow suggests that there exists a critical value n = n*, the strong disorder region for 
< n < n* and the random fixed point for < n < 1. These differences are largely 
controlled by the behavior of the structure constant Cf £ which encodes the selection 
rules in the pure 0(n) model. We interpreted the OPE intuitively and explained the 
picture of the strongly-coupled phase. Then we perform the two loop calculation near 
n — 1, where C^ £ is 0(e). The beta function is equivalent to that of the Gross-Neveu 
model and thus checked the previously known result in the random-bond Ising model. 
The correction to scaling dimensions of the energy operator and the spin operator was 
calculated up to two loop. Finally, We calculated the effective central charge and see 
this increase under the RG flow. 
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Appendix A. The vertex operators in the critical Liouville field theory 

The correlation functions in the 0(n) model can be described by the vertex operators 
in the critical Liouville theory f. The basic idea is to deform the gaussian free field 
theory by putting the background charge (— 2a ) at the infinity. This is accomplished 
by coupling the scalar field (p to the scalar curvature R. The action of the theory is 
formally given by 

S=-^- [ d 2 x UVif) 2 + 4ia R<p + \ + e ia+ * + X.e 1 ^] . (A.l) 

l07T J 

The corresponding central charge is determined from the OPE of the stress-energy tensor 
and turns out to be 

c = 1 - 24ajj. (A.2) 

The correlation function of the vertex operators V ak = e iak<f in this theory is non-zero 
only if the overall charge neutrality J2k a k = 2«o is satisfied. The screening vertex 
operators V± = e ta±ip are incorporated into the correlation functions in order to ensure 
the neutrality. The charges of the screening operators are chosen such that they are 
marginal, because otherwise the conformal invariance of the gaussian theory would be 
violated. This condition determine the charges: 

a± — ao ± \J Oq + \. (A. 3) 

The charge a_ in the 0{n) model and in the m-th minimal model are given by 

n = — 2 cos (ir/a 2 ) and a 2 = , (A. 4) 

v ' ' m+1 

respectively. In this paper, we fix our convention concerning the sign of the charges such 

that ao > and a® < describe the dilute (critical) and the dense (low-temperature) 

phase of 0(n) model, respectively. The charge a TyS and its conjugate is defined as 

a linear combination of a + and a_: 

1. 1. 
a r ,s = ^l 1 - r ) a + + ~ s ) a ~' 

a rj = ^(1 + r)a + + ^(1 + s)a_ = 2a - a TyS . (A. 5) 

The scaling dimensions of the primary operators rjS is given by 

2A rjS = -2a rtS aj^ = - [(ra+ + sa_) 2 - (a+ + a_) 2 ] . (A. 6) 

In the two-loop calculation in Section 4, we use the vertex operator representations 
of the four-point functions (£(0)£(l)£(z)£(po)) for (46) and (a(0)£(l)£(z)a(oo)) for 
(67) and (68). In order to satisfy the neutrality, we should include two of the screening 
operator V- in the vertex correlation functions. The normalization factors in the four- 
point functions are determined in a decoupling limit, such that the structure constant 

I Although it is often referred as the "Coulomb gas representation" [17, 20], we call it the "vertex 
operator representation" in the critical Liouville theory [15] in distinction with the "Coulomb gas 
method" [18]. 
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C- £ is fixed to unity. The value can be calculated by the complex Selberg integral 
as well as by the method described in Appendix D; the results are the same for both 
types of the four-point functions: 

M- l = T -^ r = ^ + 0{e), (A.7) 
where w = V Q) 2 / (2 3 V 2 ) = 2.62205755 • • • is the lemniscate constant. 



Appendix B. Some details of the integral for the two-loop beta function 
and the scaling dimension of the energy operator 

In this Appendix, we thoroughly use the result of Appendix D. We read off the exponents 
in (D.l) from the vertex operator representation (46). 2a = 4ai 3 a:Y 3 -, 26 = 8ai3aYs, 
2a' = 2b' = 4a 13 Q;_, 2f = 4a7 3 -a:_,2g = 4a 2 For convenience, we interchange the 
values of a and b, and those of a' and b' (the latter is trivial, since here a' = b'). These 
interchanges are possible by an obvious symmetry between and 1 in the integral (D.l). 

As a result, we have a = 2b = —2f = Aa^a^ = — 2 + 8e, a' — b' — —g = 2a 13 a„ = 
— §+2e. From these parameters, we determine the integral by the "scattering amplitude" 
formula (D.7). The scattering matrix M. are given by (D.8) and read, 

/ -4vre 4vre \ 



M 



+ 0(e 3 ). (B.l) 



4vre - 16vr 2 e 2 -4vre 

y -Aire + 167r 2 e 2 Ane - lQn 2 e 2 -Aire J 

From (D.20), the initial and the final state basis in (D.7) are given by 



Jt 



Jt 



4 = ( \ + 2e 2 - 4e § - 2e ) , J 3 ~ = ( \ - 6e f - 2e 1 - 4e > . (B.2) 



where, (I = 1,2,3) are real triple integrals defined in (D.9)-(D.14), and the symbol 
(...) is explained in (D.19). We should discuss the order of these integrals in e; 
for this purpose, we can use the series representation derived in Appendix D.3. From 
(D.22), the prefactors 7^ are calculated as, 

(7i + , it, 7 3 + ) =("!£ + 0(e% -16* + 0(e), -f e + 0(e°) ) , 

(B 3) 

(7i" 7a" 7 3 ") = ( - f e + 1J t + 0(e),87r + 0(e), g + ^ (-1 + 8 log2) + O(e) ). 

Naively, we expect the leading terms in J; = 7^ is the same order in e as the 7;, 
since all the triple series Si contains the (i,j,k) = (0,0,0), which is unity and hence 
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0(1). Actually, it is not the case for J+; the series S± is 0(e) because of the non-trivial 
cancellation between the constant terms. 

Since this observation makes the crucial point in the calculation, we describe the 
calculation of in detail. First, from (D.23), we have 

o+ _ y y y (-l+4eMl-4e)j(f-2e) fc (| + 2e) 3 - +fc (^ + 6e) <+J - +fc (-l + 8e) i+J - 
^^fb l] kl (l+4e) j+fc (2+4e) 

i+j+k 

Denoting the summand in (B.4) as Sij t k, we observe, from the definition of the 
Pochhammer symbol (x)k = x(x + 1) • • • (x + k — 1), that Sij^ gets a factor of 0(e) 
whenever each of two conditions {i > 2, i + j > 2} is satisfied. Thus, we let 

q i — ^0,0^) q 2 = 'S'l.o.fc) q 3 = y^ 'So,i,fc? 

fc>0 fe>0 fe>0 

«4 = ^ S'oj^, a 5 = ^ s i,j,k, (B.5) 
i>2,fc>o j>i,fe>o 

where the leading order of and 04,05 are O(e ) and (^(e 1 ), respectively. To be 

specific, at the accuracy of 0(e), we are left with the following: 



l + 4e 2 + 4e ' / 12vr 3 



4 2 y d Z V l + 4e 3 + 4e7 4 3tt 3 

a3 = ( 5 _ « ) , Fl (i - 2 « I + * \ + &\ = £ . inil! + 0(£)> (B . 8) 

3 V 8 4 y 3 V 2 + 4e 3 + 4e/ 8 9tt 3 v ; ' K J 

where the values of generalized hypergeometric functions 3F2 at unity are used. The 
arguments, which is always taken at unity, are suppressed in the first equalities in (B.7), 
(B.8) and henceforth. Now, we see the aforementioned cancellation at O(e ), and get 

cii + a 2 + a 3 = (9(e). (B.9) 

This is an example of identities known as the "contiguity relation" between hyper- 
geometric functions. Note that this type of the leading order cancellation occurs only 
for Si . We note that the cancellations at this order guarantees a consistency in the RG 
scheme. Resulting 0(e) term can be evaluated numerically, using 

(x + e) k = (x) k • {1 + e [i){x + k) — ^(x)} + 0{e 2 )} , (B.10) 

where ip(x) is the di-gamma function §. The other terms contributing (9(e) in are 

« + <* = TElt (l),(3W4), + ,(3), + , +0(£) - (B - U) 



§ Actually, there have recently been extensive studies on the closed form evaluation of the expansion 
of the hypergeometric series in algorithmic approach, for example [39]. However, to our knowledge, the 
desired expansions here have not been published. 
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To evaluate this double series, it is useful to note the following identity: 



(a)i(b)j(c) i+j (d) 



E \ a H\ u 
. flUl 



i+j 



i,3 



f (lMlMeW/) 



3-^21 



i+j 



a + b c d 
e f 



)■ 



(B.12) 



This holds because (a + b) k has the same binomial expansion as (a + b) k . By comparing 
the order 5 terms of the following expression, 

a). <-i+*>, oa« (a- 



E 



'i+j 



± + 5 =± i 



3^2 I 



(l) i (l) i (2) i+i (l) i+i 

one can express the right hand side of (B.ll) in terms of a single series: 
a 4 + a 5 = Q 1 -e + C(e 2 ), 



(B.13) 



Qi = zF, 



13 3 
2 2 2 

3 2 



-i i 

2 2 
2 1 



(B.14) 



where Fg(...) denotes the coefficients of the order 5 terms in the 3 F 2 function at unity. 
In this case, 



3-T2 



-1 1 
2 2 

2 1 



— 3-^2 1 



i + s 



-1 1 



2 1 " 2 2 
2 1 



+ C>(5 2 ). (B.15) 



We use this notation hereafter. Now, collecting the 0(e) parts of (B.6)-(B.8) and (B.14), 
we obtain 

-3tt 2 



Jt = it ■ St = + 0(e )) (0 • e° + Ae + G(e 2 )) = A + 0(e), 



(B.16) 



where we have used (B.3) for 7^. The constant A is given by 



-3tt s 
16 



15 



3^2 



3 3 1 
2 2 2 

1 3 



49 



4 

1 + 65 



3F2 



3 5 1 
2 2 2 

2 3 



± + 65 



3 3 
222 

1+45 2 + 45 



3 3 
2 2 2 

1+45 3 + 4(5 



§ -2(5 f + 25 | + 65 
2 + 45 3 + 45 



+ Qi 



(B.17) 



This completes the calculation of J+ . Since the linear relations (D.24) are solved as 



1 



J+ = - (Jf - 2Jf + J 3 -) - 2vre (jf + J 3 ") + 0(e), 



(B.18) 

J+ = — J- - 47re (Jf + J 2 -) + (9(e), (B.19) 
J 3 + = Jf _ J" + 4 7re (- J 2 " + J 3 -) + 0(e), , (B.20) 

the knowledge of the other two bases are sufficient in the formula (D.7). Actually, we 
can calculate Jf and J 3 ~ in similar manners. Let us define constants D and F as follows: 

(-1 + 4e) i (§ - 2e) i (§ - 2e) fc (§ - 10e) i+Jfc (± - 6e) i+j+fe (-| + 2e) l+J 



E 

i,j,k=0 

UJ + De + 0(e 2 ) 



fc! 



(| ~ 6e) j+fc (§ - 2e) i+i+fc (2) 



i+j 



8tt 2 

7i _ ■ 



r(l) 4 

167TC 



+ J D + C(e), 



(B.21) 
(B.22) 
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s - = (-1 + 4e) i (l - 4e) J -(| - 2e) fc (f - 10e) J+fc (| - 6e) i+J - +fc (-4e) i+J - 

3 ^ <! J' fc! (2 - 8e) i+fc (3 - 8e) i+j+k (2 - 8e) i+j 

= -^r + F€ + (<?)> (B.23) 

^ = 7 3 " ' 5 S " = ^ + F + 0(e), (B.24) 
where, again, (B.3) has been used for 7-f and 7^. Then, D and F are given by 

OF (I) 4 _2 / I 3 3 x o 2 



64 V 2 + 45 3 / 8 V 1 + 45 2 



where <2 2 denotes well-converging double series 

a -2 rf+J + , ( 3) j( i W i), (i)^ ' (B ' 27) 

Now, by substituting (B.22), (B.24) and the matrix elements (B.l) into the formula 
(D.7), we obtain 

(-2A/T 1 X(oo, e) = ^|^(2A + L> + F - J+ - J 2 ) + 16vtV ^^j^ + 

r(i) 4 

= -±±L(2A + D + F) + 0(e), (B.28) 

where, in the first and the second line, we have used (B.19) and (B.20), respectively. 
Our concern is to evaluate this, from (B.17), (B.25) and (B.26), within high numerical 
accuracy. To this end, we have used asymptotic behavior of the di-gamma function for 
the evaluation of F$. We obtained the value 

2A + D + F = 27.50074327(21), (B.29) 

which is in good agreement with another value r(l/4) 4 /27r = 4w 2 = 27.50074327208 • • •. 
Assuming the latter value in (B.28) and using (A. 7), we obtain 

J(oo,e) = -2A/'^^ + 0(e) 

8vT 

= -27r + C(e). (B.30) 



2 



Appendix C. Integrals for the correction coefficient of the spin field 
dimension 

We here calculate the integrals ^(r, e) in (64) and K,$(r, e) in (66). From (67) and 
(68), we read off the exponents defined in (D.l). After interchanging the values of the 



Disordered 0{n) Loop Model and Coupled Conformal Field Theories 



38 



V2 



parameters a <-> b and a' <-> 6' (as in Appendix B), we have a = 2a 2 3 — 2A g = | + 2e, 
6 = 2a^3T^a:i3 = — \ + e, a' = / = 2«i 3 a_ = — § + 2e, 6' = 2a^^a- = \ — e, and 
g = 2a 2 = | — 2e for /C2; we have almost the same set of the parameters as /C 2 but 
a = 2a 2 3 + 1 — 4A £ = \ + 6e for /C 3 . We substitute these into (D.8) and get, both for 
/C 2 and /C 3 , the same matrix: 

/ 2ne -27re 2vre \ 

-l + 37re 2ne -2ne + £>(e 2 ). (C.l) 
^ -1 + 3vre -1 + 37re 2vre y 

It is useful to adopt a temporary notation [x] which takes the value x and x + 4 for /C 2 
and /C3, respectively. From (D.20), the basis in (D.7) is then determined as 

/ -1 + [8]e f - 6 § - 2c v / § - [7]e f + 6 § - 2« 

J+ = ( -| + [6]e §-2e ±-e ) , Jf = ( 2 - 4e -± + 2e -I + e ) , (C.2) 



■# = ( "l + [6]e -| + 2e -± + e >,J 2 - = ( \ - [6]e -| + 2e -± + e ) ,(C3) 



-1 + [8]e 


5 
4 


— e 




-I + [6]c 


5 
2 


-2e 




1 + [2]6 


1 
2 


+ 2e 




-1 + [8]e 


5 
4 


— e 


-| + 26 


-1 + W 


1 
2 


+ 2e 




-| + 26 


1 
2 


+ 2e 


4--C 


-1 + [8]e 


3 
4 


+ e 


|-26 , 


± + 2e 


1 
2 


+ 2e 


4 + ^ 


-i + 2e 


5 

2 " 


-2e 


4 + - 



3 
4 


[7]6 


3 
4 


+ e 


|-26 


2 - 


4e 


1 
2 


+ 2e 


4 + ^ 


1 

4 ~ 


- e 


5 
2 


-2e 




3 
4 


[7]6 


5 
4 


— e 


-|+26 


1 

2 


[6]6 


1 
2 


+ 2e 




1 

4 ■ 


- e 


1 
2 


+ 2e 




3 
4 


[7]6 


5 
4 


— e 


1-26 


1 

2 


[6]6 


5 
2 


-2e 




7 
4 


[7]e 


1 
2 


+ 2e 





Jt={ | + 2e -i + 26 -i + 6 ), J 3 - = ( |-[6]6 |-26 |-6 MC.4) 

\ -1 + 26 |-26 -I + c/ \ J - [7]e -i + 26 -J + e/ 

For convenience of the presentation, we use the following notation both for /C 2 and /C 3 : 

(x, y, z) = ( j+ j+ j 3 +), (c/, v, W) = ( jf , j 2 -, j 3 -). (c.5) 

From (C.2), (D.22) and (D.23), we first observe that X is C(e°); this is guaranteed by 
the same contiguity relation as in (B.9). As a result, we realize that all of X, Y and Z 
are regular in e. By the combination of this observation and the linear relations (D.24), 
we infer V = 0(e), which suggests a non-trivial cancellation at (9(e _1 ) occurs in the 
triple series S7 determined from (C.3) and (D.23). Now, by noting V = 0(e), the linear 
relations (D.24) can be casted as 

X-Z = W/V2 + 0(e), (C.6) 
Y + Z = V2ire(U + W) + C(e 2 ), (C.7) 
fivre Y = V2ne(U + W) -V/V2 + 0(e 2 ), (C.8) 

where we have introduced the variable Q, which takes the values: 

f 8 for JC 2 

n = < (c.9) 

[12 for/C 3 . 

It should be noted that the leading order part of both U and W are 0(1), and are 
common for /C 2 and /C 3 . Using the formula (D.7) for (C.l) and (C.5), with the help of 
(C.6)-(C8), we obtain 

/C, = 27ie(U 2 + W 2 -QYZ) +C(e 2 ), (CIO) 
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where the notation Q in (C.9) is used to express both /C 2 and /C 3 in parallel. Actually, 
only the combination U 2 + W 2 is necessary for the disordered model (M — > 0). 
Substituting the parameters in (C.2) and (C.4) into (D.22) and (D.23), we obtain 

(tyi+j+k {4)1+3 , n x 

w = 97r?/2 (|)i(j)j(-j)fc (|)j+fc (|)»+j+fc rri2 x 

We obtain numerically 

f/ 2 + ^ 2 = 671.0 ±0.3, (C.13) 
which is compatible, within error bar (±0.05%), with the value 

TVIU2 6 

-Al^ = 64 = 670 .78 ... . (C.14) 

87T 6 7T 3 V ' 

Unfortunately, the convergence of the triple series U is slow, and thus our numerical 
accuracy is not good here. Nevertheless, we assume, for U 2 + W 2 , the value in (C.14). 

Appendix D. Integrals 

In our calculation of the RG functions, we should deal with a multiple integral over C 3 : 

d 2 zd 2 ud 2 v \z\ 2a \l - z\ 2h \v - z\ 2f \u - z\ 2f 

2a '\l-v\ 2b '. (D.l) 

Since this form of the integral comes from the correlation functions of the vertex 
operators (or, more plainly, from the interaction between charged particles in a two- 
dimensional plane), it seems to be ubiquitous in physics and mathematics. For example, 
the integral can be interpreted as a six-particle closed string amplitude [35] . The integral 
in a special, symmetric case of parameters (a = a', b = b' and / = g) is well studied 
in the context of twisted cohomology, and known as the "complex Selberg integral" 
[36]. What we need, however, is the formula in a non-symmetric case, when doing 
perturbation theory around a conformal fixed point. In this respect, a formula for two 
variables was used in the study of random-bond Potts model by Dotsenko et al. [20]. 
We extend their results and derive a formula for (D.l) in a systematic way. The formula, 
obtained in (D.7), takes form of a scattering amplitude. 




M 2a |i 


-z\ 2h \ 


v - 


- z \ 2 ^ \u — 


z\ 


• \u\ 2a ' 


1 1 — u 


26' 


\u-v\ 29 


v\' 



Appendix D.l. Regularization of the one- dimensional intervals 

We encounter with strong algebraic singularities that make the multiple integrals 
divergent. For this reason, we consider an analytic continuation in the parameter of 
the integrals. In order to keep the discussion clear, it is helpful to make the way of 
the analytic continuation explicit. The analytic continuation is achieved by the use of 
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J A 



a "regularization" of the intervals [36, 37], which we now describe using the following 
one-dimensional simple example. 

Consider an integral on a real interval 

-B 

dx (x- A) P (B -x) q f(x), (D.2) 

I A 

where f(x) is an analytic function on a neighborhood of the interval [A, B\. If f(A) ^ 
and f(B) ^ 0, the condition Re p < — 1 or Re q < —1 makes the integral divergent. 
The regularization of the interval [A, B] in the integral (D.2) is given by a replacement 

reg: [A, B] i — ► + [A + 5, B - 5] + (D.3) 

1 — exp{2mp) 1 — exp(2mq) 

Here, as in figure Dl, 8a and 5b are positively oriented circles of radius 5 which have 
centers A and B, and start at A + 5 and B — 5, respectively. By replacing the interval 
[A, B] by the regularized one "reg [A, B] n and taking the limit 5^0, the value of the 
integral (D.2) remains same for Rep > —1 and Reg > — 1 since the contributions from 
two additional circles vanish, and is now finite also for Rep < —1 or Reg < —1 unless 
p G Z or g G Z. In the latter case, adding two circles corresponds to the subtraction 
of infinite quantities, and the resulting finite value is what is known as the "Hadamard 
finite part" of the integral ||. 




Figure Dl. (a) Two circles 5a and 5b in the definition (D.3) of the regularization 
of a interval [A, B] . (b) The same regularization but in a simpler representation 
used in Appendix D.4. 



The definition here is natural in the following sense. Starting with the multiple 
integral (D.l) and introducing an infinitesimal imaginary part, we shall reach an iterated 
integral in which each integral is on the regularization of the interval rather than on 
the usual real interval. This regularization comes from the pairing of the paths on 
the upper and the lower half planes both of which detour the branch points. All the 
one-dimensional integrals in this Appendix should be regarded as the integral over the 
regularization of the interval. 

Appendix D.2. Decomposition of the multiple integral 

We shall consider the integrand in (D.l) as a function defined on C 3 x C 3 rather than on 
C 3 ; writing z — Z\ +iz 2 in the variable z, to be specific, the first factor in the integrand 

|| Actually, the definition (D.3) of the regularization of a interval [A, B] is proportional to well- 
known "Pochhammer contour" which is used to define the analytic continuation of the hypergeometric 
functions. 
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in (D.l): 



1 2a 



[z\ + z^Y is now denned on (z±, z 2 ) G C x C. Then the path of z 2 is 



rotated by the angle 7i/2 — 2i], without hitting any singularity, for a positive infinitesimal 
number 77. By introducing a new variable zq, we rewrite \z\ 2a as follows: 

\ z \*> = ( z j + z lf [z\ + {ie-^z,) 2 ) a 

= (z+ - irj[z + - z_]) a (z_ + ir][z + - z_]) a , (D.4) 

where the notation z± = Z\ ±z is introduced. In the following, we call {z + , u + , v + } and 
{z^, U-, V-}, respectively, holomorphic and antiholomorphic variables. Using a notation 
X v = r][X + — X_\ for a quantity X, we decompose the integral (D.l) as 

/ = ( - ) JJJJJJ dz + du + dv + dz_du-dv^ J + ■ 

[— co,oo] 6 



(D.5) 



J- 



+ 



J- 



K z + — iz ri 



■lU r 



(z^+iz v ) 



■lU r 



'(«+ 



■l+iz r . 



-u + —i[z — u 
-v + — i[u — v 
-U-+i[z — w], 



-%\z — v\ 



v y(v + -w v ) a (v+- 

^f(z~ — V-+i[z— v] 



l-iv v ) 



b' 



(D.6) 



• (m_ +i-u, ) ) a ' (u_ - 1 +iu v ) b ' (u- —v- +i [u -v] v ) g (v^ +iv v ) a ' (v- - 1 +iv v ) b ' . 

Here, J + and J71 are weakly dependent, through an infinitesimal number 77, on the 
antiholomorphic (— ) variables and the holomorphic (+) variables, respectively. If we 
fix the holomorphic variables first, the dependence of J- on the holomorphic variables 
determines relative positions of the integration paths and the two branch points and 
1 on the complex planes of the antiholomorphic variables. For this relative positions, 
we can observe there are two very distinct cases. 



(a) 




— - D - — ' E 




Figure D2. Deformation of the paths on the antiholomorphic planes. Both figures 
correspond to the sequence of the holomorphic variables z + < v + < u + . (a) Before the 
deformations. (b) After the deformations. The five symbols A, C, D, E and F 

are assigned to each regularization of the interval. The case with w_ < u_ is omitted 
here, and thus the regularization B in Figure D3 does not appear. 



The first case is trivial, and corresponds to a choice of the holomorphic variables 
such that at least one of the z + , u + and v + lies outside the interval [0, 1]. Then, according 
to (D.6), one can find at least one integration path on the antiholomorphic plane such 
that all of the points 0, 1 and the other paths projected are located on the same side of 
the plane. Therefore, the path can be contracted to a point by adding a semicircle of 
infinite radius and the corresponding integral should vanish. 
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There is, however, the second case in which all of the z + ,u + and v + belong to the 
interval (0, 1). On the antiholomorphic planes, this means all the paths of and 
V- intersect the segment [0, 1] in the same sequence as z + ,u + and v + lie on the interval 
(0,1). Take a sequence: z + < v + < u + as in Figure D2-(a), for instance. Now, we 
deform the paths so that the resulting paths encircle the point 1, the real axis from the 
upper and the lower side (Figure D2-(b)). Since the path nearest to the point 1 is that 
of m_, we first deform it, then that of i>_, and finally that of z_. Note, in deforming the 
path of, say t>_, the fixed variable w_ is projected on the t>_-plane as a branch point. 

As a result of the deformation, we have a pairing of the paths in the upper and 
the lower half planes for each variable Z-,U- and V-. Each paring of the path can be 
decomposed into a sum of the integral over the regularizations of intervals (see Appendix 
D.l for the definition of the regularizations and ref. [36] for the summation of them). 

Taking the presence of the branch cuts into account, we now attach the appropriate 
factors on these regularizations. Since each factor comes as a difference of two phase 
factors on the upper and the lower half plane, it takes the form of a sin-function. In the 
following, we assign the numbers {1,2,3} for the sequences of variables {(z± < v± < 
u ±), (v± < z± < u±), (v± < u± < z±)}, respectively. 

Consequently, from (D.6), we obtain a scattering-type formula 



7=(-2).(j+ 4 J+) 



M 



4 
\Js J 



(D.7) 



with 



M 



( s(b)s(b')[s(b')+s(g+b')} 

s(b)s(b')[s(f+b')+s(f+g+b')} 



s(b+2f)s(b')[s(b')+s(g+b')] 
s(b+f)s(b'Mb')+s(g+b')] 



,(D.8) 



s(b+f)s(b')[s(b')+s(g+b')] 
s(b)s(b') 2 +s(b+f)s(f+g+b') 
ys(b)s(f+b')[s(f+b')+s(f+g+b')] s(b')s(f+b')+s(b)s(b')s(f+g+b') s(b)s(b')[s(b')+s(g+b')] 

where a notation s(x) = sin(7r:r) is used ^f. The dimension of the basis is (3!)/2 = 3, 
where 2 comes from the symmetry between u and v. The factor 2 in (D.7) is necessary 
because of the same symmetry in J + basis. Each matrix element of M. is a sum of two 
term; each term is a product of three sin-functions attached onto the regularization of 
the intervals on and t>_-planes. Further, we have defined the initial and the final 

state basis as 



pi pu 

Ji — du Id 
Jo Jo 



dv / dzz a {\ - z)\v - z) f (u - z) } H + (u,v), 



4 = 



p\ pu PZ 

du dz dvz a (1- z)\z-v) f {u- z) f H+( u,v), 
Jo Jo Jo 



(D.9) 
(D.10) 



pi pz pu 

J+ = dz du / dvz a (l-z) b (z-v) f (z-u) f H+(u,v), (D.ll) 
Jo Jo Jo 

It would be interesting if a possible relation between the scattering matrix M. in (D.7) and A4± in 
(D.24) in this paper, and the intersection matrix and the monodromy invariant hermitian form in [37] 
was elaborated. 
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Figure D3. Pairings of the paths (the regularizations of the intervals) that contribute 
the matrix elements of A4. Each branch correspond to the rcgularization (see also 
Figure D2 for the branches A, C, D, E and F). The branches B and C, for instance, 
correspond to the regularization in the u_-plane with the sequence u_ < u_ and 
U- < V-, respectively. Each regularization has a certain factor on it due to the 
presence of the branch cuts. 



and 



poo poo poo 

J- = dz dv duz a ( 

Jl Jz Jv 



"OO poo 



"OO POO 



l) b (v- z) f (u- z) f H_(u, 


v), 


(D.12) 


l) b (z-v) f (u- z) f H_(u, 


v), 


(D.13) 


l) b {z-v) f {z-u) f H_{u, 


v), 


(D.14) 



where we have used the notation 

H+(u,v) = u a '(l - u) b \u - v) 9 v a \l - v) b \ (D.15) 

H_(u, v) = u a \u - lf{u - v) 9 v a \v - if. (D.16) 

The matrix elements of M. are conveniently understood if we draw the tree-diagrams 
which show the sequences of the variables under the process of the deformation (Figure 
D3). 

Appendix D.3. Representation of the basis through triple hypergeometric series 

By a suitable change of variables, the integral in (D.9) is transformed into an integral 
over the unit cube: 

J+= fdu [dv [dz u 2+a+2a ' +2f+9 (l-u) b 'v 1+a+a ' +f (l-v) 9 z a (l-z) f 
Jo Jo Jo 

■(l-vz) f (l-uvz) b (l-uv) b ' . (D.17) 
In general, all of Jf {I = 1, 2, 3) in (D.9)-(D.14) can be cast into the same form, namely, 

j±= [du [dv [dz u x -\i-u) y ~ l v ii - 1 {i-vy'- l ~z v -\i-~zy'- 1 

Jo Jo Jo 

■ (1 - vz)~ p (l - uv~z)- q {l - uv)-\ (D.18) 
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with each different set of exponents {A, /i, z/, A', p', v',p, q, r}. In this paper, the value of 
the triple integral in the right hand side of (D.18) is denoted as a compact symbol: 

/A A' p \ 

(p vl q j. (D.19) 
\ v v' r I 

Then each of the base jf (I = 1,2,3) in (D.9)-(D.14) looks like, 






-£-6-26' 


1 + 6 


-f 


1 


-£+a+6 


1 + / 


-b' 


-1 


-a'-b'-f-g 


1+9 


-b' 




-£-6-26' 


1+6' 


-9 


1 


- £ + a' + b' 


1 + / 


-b' 


-1 


-a'-b'-f-g 


1 + / 


-6 




-£-6-26' 


1+6' 


-/ 


1 


- £ + a' + b' 


1+9 


-6 


-1 


-a -6-2/ 


l + f 


-b' 



J 2 = 



J 3 - = 

\ -1 -a -6-2/ l + f -b' I 

(D.20) 

where, for brevity, a notation £ = 3 + a + 2a! + f + 2g is used. 

Performing the binomial expansion for the last three factors in (D.18), we get the 
following triple hypergeometric series: 

J? = it ■ (D.21) 

± _ r(A)r(AQ rqw) £(^rV) 

7i r(A+A') r(/i+/i') r(z/+z/') ' 1 ' ' 

g± _ y y y (j>)M)A r )k (A)j+ fc (//)i+j+fc m . 

' <! j! k\ (\ + \>) j+k (n + n') i+j+k (v + v>) i+j > [ ■ > 

where the Pochhammer symbol (x)k = x(x + 1) • • • (x + k — 1) is used. The parameters 
{A, n, v, A', //, z/', p, g, r} is related to the exponents {a, 6, a', 6', /, g} in (D.l) as indicated 
in (D.20). This series representation is particularly useful when some of {p, q, r, A, p., u} 
are non-positive integer plus 0(e). In that case, a separation of the order in e occurs. 



Appendix D.4- Linear relations between the basis 

For generic values of the parameters {a,b,a',b',g, f}, there exist three independent 
linear relations between the triple integral basis { J + , J~}. Although in principle we 
can evaluate the coefficients in the epsilon expansion of these integrals using the series 
expressions, but in practice some of the base are happened to be difficult to expand in 
e, while the others to be more straightforward. Hence, these relations are necessary in 
our epsilon expansion calculation of the RG functions. 
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The explicit form of the relations are, 



M. 



( Jt \ 




(jt\ 


4 


= -M- 


J2 









(D.24) 



with 



M. 



M. 



s(a)s(a') 
■2s{a'+9/2)c{g/2) 

s(a)s(a') 
■2s{a'+}+g/2)c{g/2) 

s(a)s(a'+f) 
\ •2a(o'+/-+ S /2)c( 9 /2) 

/ s(ct{W2f)s(a'W+g+f) 
1 ■2 S (a'46'+/4s/2)c( 9 /2) 

■2s(o'-H>'+/+g/2)c( g /2) 

s(a+6)s(a'+6'+g+/) 
\ ■2s(a'-+*'+/-4^/2)c( S /2) 



s(aH-/)s(a') 
•2 S (a'+g/2) C ( 5 /2) 

s(a') 

[ s ( a ) s (a')+s(a+/)s(a'+ 9 +/)] 

s(a)s(a') 
•2s(o'+/-+g/2)c( ff /2) 

s(cbff>f2/)s(a'+6'4s+/) 
•2s(a'+6'-fg/2)c( fl /2) 

■[s(oH-;>f2/)s(a , +b , -4^+/) 
+ s(aff>f/)s(a'+f)')] 

s(cH^/)s(a'4*'45+/) 
•2s(a'4*'+s/2)c(g/2) 



s(a+2/)s(a') 
•2 S (a'+g/2) C ( g /2) 

s(a+/)s(a') 
■2s(a'+g/2)c(g/2) 

s(o)s(o') 
•2 S (a'+g/2)c(g/2) / 

s((H*f2/)s(a'4*'-h?) 
•2s(a'+6'-h ? /2)c( ff /2) 

s((H*f2/)s(a'+f)'-H 9 +/) 
•2s(o'-+*'4s/2)c( S /2) 



(D.25) 



(D.26) 



s ( Q+(H .2/) S (a'-K)'+ g +/) 
•2s(a'4*'+/-+«/2)c(g/2) / 



where we have used the notation c(x) = cos(7rx). We now derive the first one of these 
three relations for an illustration. Since it is the analyticity of the integrand on the 
region except the branch cuts that makes the relation valid, the basic strategy is to 
deform successively each integration path defined on each complex plane. To keep track 
of the successive deformations of the integration paths on the three complex planes, we 
use a simple semicircular diagram to represent the regularization of the interval (see 
Figure Dl). 

Consider the z-plane in which all the branch cuts are taken along the real axis from 
each branch point to positive infinity. Since the integrand has no branch points except 
the non-negative real axis, the integral along the contour C in Figure D4 is zero. 





Figure D4. The contour C denned on the z complex plane. 



Then the ^-integral on a regularized interval reg [0, u] can be expressed by a certain 
linear combination of the integral on the intervals reg [u,v], reg [v, 1] and reg [1, 00] as 
shown in the first equality in Figure D5. Each coefficient comes from the pairing of the 
two segments which lie on opposite sides of the branch cuts. In this way, we can shift 
the integration path for each variables to the intervals in the positive real direction. 
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Jt- 




-e(a) 



-e(-f) 



-e(-2/-M 




s(a + 2/ + 6) 



-s(a+2f + b 



s(a)J+ +s(a + f)J+ + s(a+2f)J+] = 








-e(-2/) 
















»("+/) 






s(a+2/) 





s( a '+s) 




s(o'+s + 6 / 





s(a'+g + b'+f) 




s(a'+g + b'+f) 


X/f 




s(a' + b'+g) 
X 





s(a' + b') 




+ s(a' + b'+g+f) 
X 



■- s(a' + b'+g + f)J 2 + [s(a' + b') +s(a' + b'+g)] J 3 . 



«(<>' + &' + /) 



s(a' + 6' + / + <() 




ia' + b) 

— 



= [s(a' + b' + f) + s(a' + b' + f+g)] J 1 +s(a' + b')J 2 . 



Figure D5. The successive deformations of the paths on the three complex planes. 
The notation e(x) = cxp(7ra;) and s(x) — sin(7rx) are used. The semicircular shape 
represents the regularization of the interval (see Figure Dl). The dotted curves indicate 
the integrations over the variable z, while the bold curves indicate that of u or v. The 
vertical lines indicate the relative positions of the branch cuts induced for the other 
variables. 
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